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SUMMARY 


Analytical methods for combining flight acceleration and strain data with 
shake test mobility data to predict the effects of structural changes on flight 
vibrations and strains are presented. This integration of structural dynamic 
analysis with flight performance is referred to as analytical testing. The 
objective of this methodology is to analytically estimate the results of flight 
testing contemplated structural changes with minimum flying and change trials. 
The category of changes to the aircraft includes mass, stiffness, absorbers, 
isolators, and active suppressors. Examples of applying the analytical testing 
methodology using flight test and shake test data measured on an AH-1G heli- 
copter are included. The techniques and procedures for vibration testing and 
modal analysis are also described. 

INTRODUCTION 

Helicopter structural dynamics and rotor-induced vibratory loads impact the 
design, analysis, and evaluation of vibrations. The vibration environment can 
be generalized by relationships between critical points of airframe response and 
points of vibratory forcing. Vibration, conceived of in terms of how much these 
critical points shake, continues to be a major problem in helicopter dynamics. 
Structural fatigue, a result primarily of vibration, is of much greater conse- 
quence to the structural integrity of the helicopter and the reliability of its 
1 2 

systems. * The large number of airframe locations at which vibrations affect 
structural integrity, performance, and overall mission effectiveness complicates 
the issue. The vibration problem is also complicated by the many dynamic condi- 
tions that must be evaluated during helicopter development. The number of con- 
figuration changes in the operation of typical Army helicopters is enormous. If 
each of these changes is flight tested for baseline rotorcraft conditions repre- 
sented by gross weight, center of gravity, fuel loading, cargo loading, and 
flight conditions, then thousands of structural dynamic conditions must be eval- 
uated. Favorable vibrations at one airframe location due to structural changes 
may, in fact, degrade vibration at other points. Thus, the vibration solution 
process becomes one of engineering compromise. 



Structural dynamics analysis has not proven to be one of the most useful 
engineering tools in helicopter development.^ As conventionally practiced, 
most helicopter vibration tests provide limited information for resolving 
vibration issues. Helicopter flight vibration tests provide a direct measure 
of the actual vibration environment while airframe ground vibration tests are 
most often conducted to correlate analytical predictions of airframe resonances 
and mode shapes. If there is reasonable agreement between analysis and test, 
then confidence in the validity of the analysis is enhanced. However, if 
reasonable agreement is not obtained, then an impasse results. Vibration prob- 
lems are extremely difficult to quantify and have been solved by trial -and-error 
ground and flight vibration testing. 

The integration of structural dynamics analysis with flight performance, 
herein referred to as analytical testing, appears to offer a practical method- 
ology to the vibration solution process in helicopter development. This report 
describes analytical methods for combining flight acceleration and strain data 
with shake test mobility data to estimate the effects of contemplated changes 
on aircraft vibrations and stresses for various flight conditions and maneuvers. 
The category of changes to the aircraft includes external stores, weapons, 
cargo, changes in structure or materials of structure, absorbers, isolators, or 
active vibration suppressors. The objective of analytical testing methodology 
is to analytically estimate the results of flight testing such changes with 
minimum flying and change trials and to provide accurate and consistent dynamics 
information for reduced cost and testing time. 

The present investigation applied the analytical testing methodology in 
conjunction with full-scale helicopter ground and flight test vibration data. 

An AH-1G test vehicle was utilized in this project to provide ground vibration 
data of realistic quality. Flight test data of the AH-1G was obtained from 
another Army program on an as-available basis. The analytical testing examples 
in this report use AH-1G data to illustrate possible generic applications of 
the methodology. The authors do not suggest or imply applicability of these 
examples to the AH-1G or any other specific helicopter. The applications des- 
cribed herein were directed to the practical acquisition of helicopter vibra- 
tion data for analytical testing and the possible utilization of the method. 
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The reader is cautioned not to interpret the results of these illustrations as 
representing experimental validation of analytical testing in any possible 
appl ication. 
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LIST OF SYMBOLS 

modal acceleration, g/N (g/lb) 

change factor (a constant), such as the proportionality of skin 
thickness of the change 

generalized damping, N-s/m (lb-s/in.) 

diameter of modal circle, m/N (in. /lb) 

damping matrix, N-s/m (lb-s/in. ) 

rigid body acceleration coefficient, g/N (g/lb) 

force, N(lb), or frequency, Hz 

mode frequency function, defined in text 

structural damping coefficient or units of acceleration 

gravitational acceleration, 9.8066 m/s 2 (32.174 ft/s 2 ) 

imaginary operator, A 1 

generalized stiffness, N/m (Ib/in.) 

stiffness matrix, N/m (lb/in.) 

absorber mass or generalized mass, kg (Ib-s /in.) 

moment, N-m (lb-ft) 

mass matrix, kg (Ib-s /in.) 

rotor revolution, thus 2P denotes twice rotor revolution 
vibratory displacement, m (in.) 
vibratory acceleration, g 

vector of complex responses, defined in text 
mobility residual, g/N (g/lb) 

arc length of modal circle or, in Matrix Difference Method, 


A “x 



time, s 

transformation matrix = 

_aq 

sweep speed required for swept sine shake testing. Hz/s 

maximum speed attainable in level flight at maximum continuous 
power 
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x displacement, m (in.) 

x linear acceleration, m/s (in./s ) 

(y(w)} displacement vector of a steadily vibrating undamped multiple 
degree of freedom system 

(y(o))> Fourier transform of the accelerations 

Y displacement mobility, m/N (in. /lb) 

C 

Y mobility of structural change, m/N (in. /lb) 

Y^ strain mobility, N" 1 (lb -1 ) 

Y acceleration mobility, g/N (g/lb) 

Z impedance, N/m (Ib/in.) in Appendices A and B 
Z c impedance change 

[Z] impedance matrix, N/m (Ib/in.) 

a angle of bar in simple bar truss or logarithmic sweep speed for 
swept sine shake testing, dec/min 

6 frequency ratio 

6 distance between col inear displacements, m (in.) 

A indicates change in variable following the symbol 

A.p mode resolution, Hz 

e strain, nondimens ional 

C viscous damping factor 

9 rotation response, deg 

p 

a(o)) frequency dependent complex eigenvalue, Hz 

y frequency function 

t characteristic time for secular response, s 
<j> phase angle 

$ normalized mode 

v orthonormal mode, kg”^ 2 (in.^ 2 /lb^ 2 -s) 

( £ \ 

v' orthonormal strain mode 

w forcing frequency, rad/s 

natural frequency of the nth mode, rad/s 
Qj absorber tuning frequency , rad/s 

( )* imaginary component of complex variable 

D 

( ) real component of complex variable 
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/ \* 

( ) complex conjugate 

( )/<j> denotes amplitude and phase of complex variable 
[ ] rectangular or square matrix 
{ } column matrix 

|_ ] row matrix 

[ ] T matrix transpose 

[ matrix inverse 

[ ] + matrix pseudoinverse 

C' ^ diagonal matrix 

[“I] identity matrix 

9/9x partial derivative 




THE PRACTICAL ASPECTS OF ANALYTICAL TESTING 


Analytical testing can be applied in three stages of helicopter develop- 
ment; namely, during preliminary design, after construction but before flight, 
and during development flight testing. In other words, the vibration data can 
originate from analysis or test or in combination. 

During the first stage of helicopter development, the helicopter exists on 
paper and flight vibrations are simulated using available theoretical informa- 
tion. Hub excitation estimates are obtained from rotor analyses and estimates 
of downwash impingement and wake interferences can be obtained from aerodynamic 
considerations. These estimated excitations for various maneuvers can be 
applied to mobilities obtained from airframe finite element or similar dynamic 
analyses to yield estimates of in-flight stresses and vibrations. The predic- 
tion of helicopter vibrations ultimately requires consideration of the response 
of the coupled rotor/airframe. In general, the degree of rotor/airframe coup- 
ling is strongly dependent on the characteristics of both structural and aero- 
dynamic interface coupling. 

Substructures and changes can be combined with the basic structure in the 
analytical testing processes which can be implemented on interactive minicom- 

o 

puters. Therefore, the finite element program for each structural module 
needs to be executed only one time to yield resonance and mode shape data. 

Modal damping can be accounted for during analytical testing. This process 
minimizes finite element program sizes and running times while allowing a large 
number of nodes in the combined structure. In addition, it allows the repre- 
sentation of nonproportional damping and, consequently, complex modes. The 
complex modes, resonances, and modal damping of the combined structure can be 
obtained by treating the combined mobilities as test data. Uncertainties of 
theoretical resonances and modal damping can be evaluated by direct perturba- 
tion of these parameters in the modal synthesis. The modal synthesis may also 
include high frequency residual mobilities for local mode driving-point effects. 
Uncertainties in rotor loads estimates^ and in estimates of downwash impingement 
and wake interference can be systematically assessed with variation in excita- 
tion magnitudes and phasings during this stage of helicopter development. 
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In the second stage of helicopter development, after construction but 

before flight, a shake test aircraft should be available. This non-flying 

shake test aircraft would be used for analytical testing of changes, ground 

flying for fatigue evaluation , 5 and analysis of accidents. In this stage, 

shake testing data can be used to refine the existing finite element 
6 7 8 

models. 5 Mobilities from finite element models of substructure can be com- 
bined with mobilities from tests of other subsystems to obtain total system 
mobilities, resonances, modal damping, and complex modes. In general, the 
fuselage system mobilities and the mobilities of complex components, such as 
engines, may be obtained through shake testing while mobilities of contemplated 
changes, flexible portions of the airframe, and low mass structural appendages 
would be obtained through finite element modeling. This separation between 
finite element modeling and shake testing for analytical testing purposes 
offers optimum utilization of finite element analysis and modal analysis shake 
testing. The major accomplishment in this stage is the application of refined 
mobilities for identifying favorable structural changes to adjust airframe 
resonances and nodes. Theoretical estimates of the vibratory loads are com- 
bined with these mobilities to estimate changes in vibrations and stresses. 

During the third stage of helicopter development, analytical testing uses 
only ground and flight vibration data to estimate changes in flight responses 
caused by structural and configuration changes. As the flight envelope of the 
prototype aircraft is expanded, flight vibration data can be applied directly 
since the theoretical estimates of external excitations are no longer needed. 
Anticipated changes such as stores, weapons, cargo, and structure can be exam- 
ined for flight effects on vibrations and stresses before actual flight. Know- 
ledge of the number, types, or locations of the external excitations is not 
required. The only flight data necessary are accelerations and strains mea- 
sured during the initial flight tests. 
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ANALYTICAL TESTING THEORY 


The matrix equations of motion for a damped linear structure can be gener- 
alized in the frequency domain as 


[K] - a) [M] + i [D(w)] 


(q> = (f } 


0 ) 


where [K], [M], and [D(w)] are Nth ordered stiffness, mass, and damping matrices, 
respectively. The matrix terms on the lefthand side of equation (1) define the 
displacement impedance matrix, [Z], or 


[Z] (q> = (f > ( 2 ) 

The dynamic responses can thus be characterized by the simple matrix equation 
given by 


(q) = [Y] if) (3) 

where 

[Y] = [Z]’ 1 (4) 

and the variables in equation (3) are complex valued and frequency dependent. 

The matrix [Y] is the transfer function which relates the input excitations, 

{ f } , to the output responses, { q } . If the response vector is displacement, 
velocity, or acceleration, then the transfer function is defined as the displace- 
ment, velocity, or acceleration mobility, respectively. Compliance, mobility, 
and inertance are sometimes used in the literature to define the corresponding 
displacement, velocity, and acceleration mobilities. 

Types of Mobilities 

Mobilities are defined as partial derivatives of response with respect to 
excitation in the frequency domain. In general, there are two separate types 
of mobilities used in analytical testing which can be distinguished by consid- 
ering the nature of the response and the excitation. 
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If the response is vibration (displacement, velocity, acceleration) and the 
excitation is force or moment, then the mobility matrix consists of the fol- 
lowing components 


■[V] = 


' 3 V 3f k 8 V 9H k 

8 V 3f k 3 V 3M k 


(5) 


where q and e are the respective translation and rotation responses; f and M 
are the force and moment excitations, respectively. The reciprocity principle 
for vibration mobilities is satisfied such that 


8 V 3 f k = 3 V 3f j 


3ej / 3 f k = 3 q k /3M j 


90j/3M k = 90 k /3M. 


(6) 


If the response is strain and the excitation is force, then the strain 
mobility matrix becomes 


’ yU) ] ■ b-/ 3f k ] 


(7) 


Let q - , represent a displacement col inear with q. and separated a distance 6 . 

vJ ' I yj yj 


from q . . Then, the strain at j is defined such that 


e j s ji ! 0 (qj+i ■ q j )/s j 


(8) 


and the strain mobility becomes 


3 e j / 3 f k E ( Vi,k - Y jk )/S j 

J 


(9) 


This type of strain mobility is useful for evaluating stiffness changes. The 
application of de/df to analytical testing is discussed in the next chapter. 
Note that the reciprocity principle is not satisfied for strain mobility; i.e.. 
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(10) 


3 £j /3 f k i* 3e k /3fj 


Analytical Testing Equations 


The solution to helicopter vibration problems consists, in part, of pre- 
dicting and confirming the flight vibration or strain effects of dynamic changes 
in the airframe, on the rotor, or at the rotor/airframe interface. Let (q > and 
{e} represent the baseline vibration and strain, respectively. Then, (q ' > and {e 1 } 
represent the change in baseline vibration and strain due to a dynamic change. 

In theory, {q 1 } and {e 1 } can result from either a change in the mobility matrix 
or a change in the vibratory loading. The basic analytical testing method con- 
siders changes in the mobility matrix which can be synthesized by discrete as 
well as multiple and distributed dynamic impedance adjustments. The category 
of impedance changes that can be accommodated includes mass, stiffness, absor- 
bers, isolators, and active suppressors. As shown in Appendix A, the matrix 
equation for determining the change in baseline vibration for a general multi- 
dimensional impedance change is 


{q ' ) = (q > 



[ Y n + Y ir 



(ID 


where Yjj is defined from the impedance change. The mobilities Y^j and Y^ 
represent transfer functions for the baseline structure at the change interface 
and do not include the effects of the impedance change. In other words, for a 
properly modeled impedance change, the effects on flight vibrations are evalu- 
ated without incorporating the change in the baseline structure. Therefore, 
only one NASTRAN or similar dynamic analysis is required to implement equa- 
tion (11). If changes in strain, as opposed to vibration, are considered, then 
equation (11) becomes 

-1 


{ e 1 } = (e ) 


f 

_ A 


v( e) 

Y ql 


Y + Y 
Y II y II 


02) 


In summary, the changed flight responses (vibration or strain) are charact- 
erized by the dynamics of the impedance change, the baseline flight responses, 
and the baseline mobility responses. The operational equations can be used with 
either theoretically derived vibration and strain data or ground and flight 
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vibration and strain measurements. In the next chapter, equations (11) and (12) 
are considered for examining mass, stiffness, and absorber changes to illustrate 
the analytical testing methodology. 

Limitations of the Method 

Basic to the analytical testing method are the assumptions that the struc- 
ture is linear and that changes to the structure do not charge the external 
loadings. These conditions are only approximated in an actual helicopter 
flight. The airframe is not a linear system, as shake tests of the AH-1G 
showed, but it appears that it can be represented as linear for most practical 
engineering purposes. 

The second assumption is a workable approximation under some conditions of 
change and not under others. It is important to note that the mobilities used 
in analytical testing must be physically realizable and consistent for any 
driving-point. Mobilities obtained using a lumped mass representation of the 
rotor at the hub violate this requirement but the practical effect of the viola- 
tion is not known. Ideally, the mobilities of the airframe would contain the 
dynamic effects of a rotating rotor in a vacuum and this might be approximated 

by coupling theoretical rotor mobilities with shake test airframe mobilities in 

the analytical testing equation, except that the partial derivatives of in-plane 
hub motions with respect to in-plane hub forces have periodic coefficients. It 
is not analytically difficult to handle this problem, and therefore remove this 
particular limitation, but the practicality of doing so has not been estab- 
lished. The effect on change estimates of airframe mobilities without a 

rotating rotor in a vacuum is discussed in Appendix B. 

The basic analytical testing equation is of the form 

{q • > = {q} - [A] {r} (13) 

where q is the vector of complex motions (vibrations or strains) on the air- 
frame measured in flight without a structural change, r is the vector of com- 
plex motions measured in flight at the coordinates of the structural change, 
q' is the vector of complex flight motions that result from the change and A is 
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a matrix function of measured airframe mobilities and mobilities of the struc- 
tural change. If the structural change has negligible effect on the structural 
dynamics of the helicopter, the A matrix is nearly null and Ar is negligible; 
resulting flight vibrations and strains are virtually unchanged. If, on the 
other hand, the structural change has a significant effect on the dynamics such 
that the absolute values of the Ar terms are much greater than the absolute 
values of the q terms, then the change will make the flight vibrations and 
strains much higher. It is seen, therefore, that analytical testing is least 
sensitive to errors in mobilities or modeling of the structural change when: 

(1) the structural change has negligible effect on flight vibrations and strains; 
and (2) the structural change results in much higher flight vibrations and 
strains. 

The practical consequence is that one does not need high precision mobili- 
ties or change modeling to filter out quite rapidly those structural changes 
which either do not significantly change flight vibrations or strains or those 
structural changes which will significantly worsen flight vibrations or strains. 
Since most changes contemplated in the life cycle of a military helicopter are 
for mission improvement, not dynamics or stress improvement, there is an 
obvious value in using approximate but not precise structural dynamics data to 
identify those changes which are likely to create serious flight structural 
problems before proceeding with flight testing the changes or with more expen- 
sive and more precise measures of analysis and test. 

APPLICATIONS OF ANALYTICAL TESTING 

The following numerical illustrations of the analytical testing methodol- 
ogy utilize ground and flight vibration data obtained from an Army AH-1G test 
vehicle. The types of dynamic changes which are considered include mass, stiff- 
ness, absorber, and active suppressors. Except for the mass change example, 
the analytical testing illustrations are hypothetical. In addition, these 
results do not suggest appl icabil i ty to the AH-1G or any specific helicopter. 
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Test Vehicle and Test Conditions 


The AH-1G is an armed helicopter that is configured with two-bladed, tee- 
tering main and tail rotors. Armaments include nose mounted and wing mounted 
weapons. A plan side view of the AH-1G is shown below. 



Mass Changes 

Any change of a structural nature in a helicopter can be described in terms 
of its attachment point mobilities which include the effects of stiffness, mass, 
and damping of the change. However, the most common changes can be described 
essentially as mass changes. Among these are external stores added to or taken 
from the aircraft as exemplified by rocket pods, bombs, missiles, guns, and 
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external fuel. There is often an extremely large number of possible combina- 
tions of external stores on helicopters because of the variety of missions. 
Cargo and transport helicopters have many variations in payload of a mass 
change nature. There are mass changes from fuel burn-off, depletion of ammu- 
nition, and firing of rockets or missiles. 

In the continuing development of a helicopter it would be expeditious to 
predict the effects of flight stresses and vibration of such changes so that 
problem areas can be anticipated, engineering judgments can be made, and cor- 
rective action prepared with minimum trial -and-error testing. Allowing for 
flight data scatter, the engineer would make such predictions at several criti- 
cal locations with analytical testing for critical classes of maneuvers. It is 
impractical to attempt precise predictions for every possible airspeed, gross 
weight, c.g. location, yaw rate, pitch rate, roll rate, power setting, air tem- 
perature, wind condition, altitude, etc., and for all the locations of interest 
on the helicopter. In a well developed helicopter analytical testing would be 
used for major changes, such as the contemplated addition of rocket pods, but 
in a helicopter in the early stages of development analytical testing would be 
applied to a wider variety of mass changes to aid in identifying possible 
problems. 


A Flight Example for Mass Changes 

The AH-1G is used to illustrate the analytical testing methodology. Con- 
sider the hypothetical situation of an AH-1G which had never flown with rocket 
pods. An addition of rocket pods weighing 181 kg (400 lb) each to the outboard 
wing station is contemplated. This represents a 9.4 percent increase in the 
gross weight of the helicopter. The predictions are made from the clean con- 
figuration for classes of maneuvers without accounting for additional drag, 
fuel burnoff, or other causes of possible changes in external aerodynamic 
loading. 

The following flight acceleration data were taken on an as-available basis 
from another project. No strain data were available. The clean configuration, 
without rocket pods, had a take-off gross weight of 3830 kg (8465 lb) and was 
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flown in ground winds of 5 to 7 knots. Outside Air Temperature (OAT) of 10°C 
(50°F) and 766 mm (30.17 in.) Hg barometric pressure. The flights with the 
rocket pods were made at 4106 kg (9075 lb) take-off gross weight with ground 
winds of 3 to 5 knots, OAT of 20°C (68°F) and 754 mm (29.68 in.) Hg barometric 
pressure. The c.g. was at FS 196.3 in both flights. Power and control settings 
were not necessarily matched in the flights. Data were analyzed for the condi- 
tion of highest peak-to-peak vibration of a set of selected points, differing 
from flight to flight, in each class of flight condition with a harmonic anal- 
ysis over five to eight rotor revolutions. Except for the unavailability of 
strains, this situation is representati ve of practical application of analy- 
tical testing. 

The mobilities of the left and right wing stores positions for vertical 
motion constitute a 2 x 2 complex matrix. These are at Butt Lines +60 and are 
identified as Z200L and Z200R. 



. 054/-8 0 
.012/14° 


.01 2/ 1 4°1 

.040/8° 1000 N 


(14) 


Analytical testing predictions were made for nine motion coordinates shown in 
the sketch below. 



15 



The matrix of mobilities for these nine coordinates relative to forcing at Z200L 
and Z200R is as follows: 



Z200L 

Z200R 


Z200 

".033/ -4° 

.052/9° 


Z90L 

.017/22° 

.024/-140 0 


Z140L 

.046/5° 

. 01 8/ - 1 76° 


Z396R 

. 046/-1 9° 

.068/11° 


Z400 

. 1 00/ -49° 

. 049/-97 0 

9 

1000 N 

Y380 

.193/144° 

. 208/-48 0 

Y440 

. 091 / - 1 6 7° 

.086/8° 


Y490 

.065/77° 

. 084/-96 0 


Y51 7 

_. 168/1 30° 

. 1 51 /-53° 



change is obtained from equa- 


L 0 181 kgj 

where G is the gravitational acceleration. The effects of this mass change, 
using equation (16), on 2P(10.8 Hz) vibrations are summarized in the following 
discussions. 

Level flights . - The effect of the pods at . 5V^, shown in Table I, is 
small. The aft tail lateral accelerations, which are sensitive to main rotor 
2P forces on the vertical fin, are higher than predicted but still low. At 
the pods result in lower vertical vibration, as shown in Table II, but the 
change is not very great. The lateral vibration at FS Y380 was predicted to 
increase and showed a greater increase in the pod flight. The decrease in 
vibration at FS Y490 and FS Y517 was significantly greater than predicted. 
Whether the difference is due to inaccuracy of the method, the 181 kg (400 lb) 
weight difference, or the 8 knot air speed difference is not known but the gen- 
eral conclusions from the prediction are reflected in the flight with the pods. 
Gross weight difference results from fuel burnoff. 



The vibration, in g units, resulting from the 
tion (11) 


= q 


Y 1r 


1 


181kg 
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TABLE I. - 2P VIBRATION IN STRAIGHT AND LEVEL FLIGHT 
AT .5 V H , g 



Flight of clean 
configuration 

Predicted for 
flight with pods 

Flight with pods 

. 

3840 kg (8465 lb) 
gross weight 

4203 kg (9265 lb) 
gross weight 

4035 kg (8895 lb) 
gross weight 


70 knots 

67 knots 

Z200 

.10/0° 

.09/0° 

.11/14° 

Z90L 

.03/83° 

.02/73° 

.03/100° 

Z140L 

.04/33° 

.03/24° 

.03/43° 

Z396R 

. 25/-21 ° 

. 23/-22° 

. 21 /-9° 

Z400 

. 22/-1 9° 

.21/-16° 

. 20/-8° 

Y380 

.07/262° 

.12/85° 

.10/119° 

Y440 

.14/233° 

.16/126° 

.21/137° 

Y490 

.18/22° 

.17/133° 

.26/148° 

Y51 7 

.29/201° 

.29/151° 

.37/167° 


TABLE II. - 2P VIBRATION IN STRAIGHT AND LEVEL FLIGHT 
AT V H , g 



FI ight of clean 
configuration 

Predicted for 
flight with pods 

Flight with pods 


3840 kg (8465 lb) 
gross weight 

4203 kg (9265 lb) 
gross weight 

4017 kg (8855 lb) 
gross weight 


144 knots 

136 knots 

Z200 

.20/73° 

.18/73° 

.15/85° 

Z90L 

.20/133° 

.18/129° 

.18/125° 

Z140L 

.18/114° 

.15/109° 

.16/104° 

Z396R 

.31/48° 

.27/48° 

.22/48° 

Z400 

.34/44° 

.32/45° 

.31/46° 

Y380 

.19/206° 

.25/156° 

.31/177° 

Y440 

. 61/-145° 

. 66/-1 52° 

. 64/ -1 67° 

Y490 

. 88/-1 50° 

. 82/ —1 52° 

. 50/-1 48° 

Y51 7 

1.66/-138 0 

1 . 57/-1 42° 

1 . 1 5/ -1 44° 


Gunnery runs . - The data for the flight of the clean aircraft were not 
necessarily taken for a portion of the rolling pullout comparable to that for 
which those data were taken on the flight with the pods. From the prediction 
of the effects of the pods, shown in Table III, the pods have little effect on 
the vibration at some locations and cause a decrease in the vibration at others 
in a high load factor rolling pullout to the left. The pod flight data leads 
to the same conclusion. The same situation pertains in a rolling pullout to 
the right, shown in Table IV. In this case the predictions for FS Z396R and 
FS Y517 showed reductions that were not seen in the pod flight as did the pre- 
diction for FS Y380 in Table III. However, considering the nature of the data 
samples, the conclusion of the prediction is reflected in the pod flight. 

Sideward flight and landing . - The flight data in Table V are samples 
taken from sideward flight to the right to 35 knots with reversal and sideward 
flight to the left to 35 knots. There is no significant change in vibration 
due to the pods in sideward flight or in approach and landing. The phase angles 
of the data at very low g-levels, determined from harmonic analysis, are highly 
variable. In approach and landing with pods. Table VI, the data sample shows 
low vibration at FS 490 and FS 517, fin stations of large vibration scatter, 
and is most likely the result of the time sample chosen. 
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TABLE III. - 2P VIBRATION IN ROLLING PULLOUT TO LEFT, g 



Flight of clean 
configuration 

Predicted for 
flight with pods 

Flight with pods 


3645 kg (8035 lb) 

4130 kg (9105 lb) 

3966 kg (8744 lb) 


gross weight 

gross weight 

gross weight 


187 knots 

186 knots 


1.42 g load factor 

1.52 g Toad factor 

Z200 

.53/85° 

.52/89° 

.58/84° 

Z90L 

.46/118° 

.42/117° 

.42/116° 

Z140L 

.47/105° 

.42/103° 

.40/104° 

Z396R 

.52/73° 

.42/75° 

.54/76° 

Z400 

.94/68° 

.83/70° 

.79/66° 

Y380 

. 50/-1 26° 

. 32/-1 58° 

. 52/ -1 57° 

Y440 

1 . 28/-11 5° 

1 . 25/-1 20° 

1 . 04/ -1 41 ° 

Y490 

1 . 27/- 1 16° 

1 . 19/-1 14° 

. 93/- 127 ° 

Y51 7 

1 .82/-91 ° 

1 .60/-89 0 

1 . 59/-92° 


TABLE IV. - 2P VIBRATION IN ROLLING PULL 

.OUT TO RIGHT, g 


Flight of clean 
configuration 

Predicted for 
flight with pods 

FI ight wi th pods 


3792 kg (8361 lb) 

4155 kg (9161 lb) 

4042 kg (8909 lb) 


gross weight 

gross weight 

gross weight 


164 to 128 knots 

162 knots 


1.78 to 1.5 

g load factor 

1 .23 g load factor 

Z200 

.36/84° 

.32/85° 

.24/79° 

Z90L 

.49/129° 

.44/128° 

.31/117° 

Z140L 

.43/115° 

.37/113° 

.31/106° 

Z396R 

.33/55° 

.26/56° 

.47/71° 

Z400 

.82/61° 

.72/61° 

.63/60° 

Y380 

. 47/- 121 ° 

.24/193° 

.24/-167° 

Y440 

. 97/- 1 1 4° 

. 94/ -1 23° 

. 79/-1 58° 

Y490 

. 71/-121 0 

.61/-1 1 6° 

.47/-1 44° 

Y 51 7 

.85/-106 0 

. 58/-1 07° 

. 85/ —1 23° 


TABLE V. - 2P VIBRATION IN SIDEWARD FLIGHT, g 



Flight of clean 
configuration 

Predicted for 
flight with pods 

Flight with pods 


3688 kg (8130 lb) 
gross weight 

4051 kg (8930 lb) 
gross weight 

3980 kg (8774 lb) 
gross weight 


35 knots right to 35 knots left 

Z200 

.23/142° 

.20/143° 

.25/130° 

Z90L 

.11/145° 

o 

m J 

•-J 

o 

.12/138° 

Z140L 

.13/144° 

.12/144° 

.12/138° 

Z396R 

.26/137° 

.22/136° 

.18/131° 

Z400 

.25/132° 

.22/140° 

.23/125° 

Y380 

. 03/-46° 

- 07 / -31 0 

.02/-52 0 

Y440 

. 08/-43° 

. 06/ -47° 

.08/-82° 

Y490 

.03/54° 

.03/47° 

.04/104° 

Y51 7 

.02/45° 

.04/73° 

- 01 / -1 42 ° 


TABLE VI. 2P VIBRATION IN APPROACH AND 

LANDING, g 


Flight of clean 
configuration 

Predicted for 
flight with pods 

FI ight with pods 


3613 kg (7966 lb) 
gross weight 

3976 kg (8766 lb) 
gross weight 

4002 kg (8824 lb) 
gross weight 

Z200 

. 1 4/-1 3° 

. 1 3/-1 3° 

. 1 4/-21 ° 

Z90L 

.10/31° 

.09/29° 

.07/24° 

Z140L 

.11/12° 

.10/9° 

.08/2° 

Z396R 

. 24/-47° 

.22/-49° 

.20/-38 0 

Z400 

. 27/-24° 

. 25/-20 0 

.25/-28° 

Y380 

.10/126° 

.10/82° 

.08/111° 

Y440 

.18/140° 

.19/131° 

.15/-121 ° 

Y490 

.15/137° 

.13/136° 

o 

1 

''J 

-p» 

0 

Y51 7 

. 24/-40 0 

— — 

. 29/-33° 

.05/ -1 61 ° 
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Vibration Absorber Changes 


All major helicopter manufacturers have used conventional spring-mass 
vibration absorbers with varying degrees of success. Except for in-plane cen- 
trifugal hub absorbers which are used to cancel (N-l )P and (N+1)P shears in the 
rotating system of helicopters, absorbers are customarily tuned to be coinci- 
dent with the excitation frequency and placed at the airframe location where 
vibration reduction is needed. In limited applications, airframe absorbers 
have been remotely placed from points where low vibration is desired. The 
remote absorber has the advantage of providing vibration reduction at locations 
where conventional airframe absorbers are physically impractical. These remote 
absorbers have been trial -and-error tuned to be somewhat off resonance to 
achieve optimum vibration reduction. 

Airspeed, gross weight, and center of gravity variations alter the rela- 
tive magnitudes and phasings of airframe responses, thus absorber effectiveness 
varies. There are two situations in which the effects of absorbers are indepen- 
dent of the airspeed and flight maneuver of the helicopter. First, is the 
extreme situation with absorbers coinci'dent to each vectorial coordinate of 
external forces and moments acting on the helicopter. In this case, the entire 
airframe has virtually zero vibration at all airspeeds and in all maneuvers. 

The second situation is at the attachment point of the absorbers along the 
absorber direction. Clearly, the dynamicist must consider the flight effects 
of a limited number of absorbers at airframe locations remote from the 
absorbers . 


For single-point or discrete impedance changes, equations (11) and (12) 
reduce to simple scalar algebra equations for estimating the effects on vibra- 
tions or strains. The vibration absorber analytical testing equations become 
for acceleration, after simplifying equation (11), 


Y. q 
jr t 

Y c + Y 
rr rr 


(17) 
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and, for strain, after simplifying equation (12) 


Y>> q 

■ = _ jr H r 

J J v c + Y 

rr rr 


( 18 ) 


where Y r £ is the unrestrained driving-point acceleration mobility of the absor- 
ber at its attachment point r. For a structurally damped absorber, the unre- 
strained driving-point acceleration mobility is 



2 2 

l-o) /fiy + ig 
m(l+ig) 


(19) 


where fty is the tuning frequency, m is the absorber mass, and g is the non- 
dimensional structural damping coefficient. Similarly, with viscous damping 

2 2 

1-0) /QZ + i2£0)/Q T 
Y rr " m(l+i2co)/Oj.) 

where ; is the viscous damping factor. As equations (17) and (18) show, the 
required absorber weight, tuning frequency (not necessarily equal to the excita- 
tion frequency), and damping (not necessarily zero) for minimum vibration or 
strain along the motion coordinate j are functions of the flight vibration at r, 
the flight vibration or strain at j, the r-j mobility, and the r-r mobility. A 
vibration absorber can be designed to minimize both vibrations and strains. For 
vibratory strains the required r-j mobility is defined as the strain mobility. 


Vibration attenuation is a fundamental objective in dynamics and equa- 
tion (17) can be written in a slightly different form as 


= 1 *° ' 


Y. q /q. 
jr H r 

Y c + Y 
rr rr 


( 21 ) 


The effects of tuning and damping on vibration at j can be examined by plotting 
the absolute value of (q'./q.) versus the tuning ratio, fiy/w, for fixed values 
of damping. Interactive computer graphics offer the capability to investigate 
the effects at various critical airframe locations to select the tuning fre- 
quency for the optimum overall effect. A weighting factor can be assigned to 
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|q'/q| for each station and one composite curve can be displayed for many 
stations. 


Conventional absorbers are analyzed with regard to vibration at the absor- 
ber attachment point r. If j replaces r in equation (21), then the absorber 
transmissibil ity becomes 





+ Y. 


rr 


( 22 ) 


or, after substituting for Y^ r in equation (19), 


q^/q r = 


? ? 

1 - oj /fi j + i g 


1 - ^ 

n 


2 + ^rr m " ^rr mg + 1 
T 


’ ^Y 1 m + mgY R + g 
rr s rr y 


(23) 


The percentage reduction of vibration at the absorber attachment point is, of 
course, independent of the flight vibration. A conventional absorber with zero 
damping produces zero vibration when Y^ r is zero or when the tuning frequency 
coincides with the excitation frequency. When the forcing frequency equals the 
tuning frequency, equation (23) gives 


q^/q 


r 


= £ 
m 


- 



1 



g + i 

Y^ + g Y^ + 

2-1 

rr rr s 

rr 3 rr 

mj 


Thus, for a given absorber damping the attachment point vibration ratio is 
inversely proportional to the absorber mass. 


(24) 


The imaginary part of the helicopter driving-point acceleration mobilities 
is necessarily positive, but the real part may be either positive or negative. 
Because helicopter structures do not have proportional damping the modes are, 
in general, complex and, even in the vicinity of a well separated mode with a 
high rr modal acceleration (residue), the real part of the driving-point mobil- 
ity is not necessarily smaller just below resonance than just above a resonance 
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as in structures with classical modes. For structures with classical modes, 
the signs of the real parts of the driving-point mobilities may be the same 
just above and just below a resonance, due to coupling of other modes, and 
that sign may be positive or negative. Therefore, the resonance introduced 
by the addition of the absorber may occur above or below the driving-point 
anti resonant frequency created by the absorber. Because the driving-point 
mobility is a function of frequency it is necessary to examine equation (22) 
over a frequency spectrum to determine the absorber bandwidth, that is, the 
change in resultant vibration with variation in excitation frequency or rotor 
speed. 


As seen from equation (24), the vibration ratio at the attachment point 
is, for small damping, directly proportional to the damping. This leads to the 
misconception that minimal absorber damping is desirable regardless of the loca- 
tion of the absorber and the point of concern on the helicopter. 


Obtaining zero vibration . - In equation (21) let q'. be zero. Then, 

J 



Y 


rr 


( 25 ) 


The absorber mobility required for zero vibration along the j motion coordinate 
is a function of the in-flight vibration of the motion coordinate j and that of 
the absorber, r. The complex q’ r /q . ratio will usually be different for each 
flight condition and zero vibration at j from an absorber at r is dependent on 
the flight condition. Since the dri ving-point imaginary acceleration mobility 
must be positive, it is necessary but not sufficient for the imaginary part of 
the right hand side of equation (25) to be positive for zero vibration. 

Through control of absorber damping zero vibration at a desired point can be 
achieved in any one flight condition only under certain circumstances. 


24 



Examples of Absorber Analysis 


The following examples are presented to illustrate the applicability of 
analytical testing for examining the effects of absorbers on AH-1G airframe 
vibrations. The absorber is a simple spring-mass device with hysteretic damping 
as discussed for equation (19). Equations (21) and (22) are used to determine 
the performance of the absorber at critical points on the airframe for several 
representative flight conditions. 

Absorber at nose . - This example considers the effects on vertical vibra- 
tion at the gunner's left ( FS Z90L) and at the horizontal stabilizer (FS Z400). 

A vertical absorber, weighing 13.61 kg (30 lb), is located at the nose (FS Z50) 
of the AH-1G as shown in the sketch below. 


The acceleration mobilities which were measured during the ground vibra- 
tion test are 



Y (Z50, Z50 ) 

- Y (Z90, Z50) • 

Y (Z400,Z50 ) 


'.100 g/1000 N (.044 g/100 lb)/10° 
-.09 g/1000 N (.040 g/100 lb)/8° ■ 

.232 g/1000 N (.103 g/100 lb)/123° 


(26) 
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The vibrations measured in flight are shown in Table VII. The driving- 
point mobility at FS Z50 is nearly the same as the transfer mobility between 
Z50 and Z90L and the flight vibrations of these coordinates are only slightly 
different in magnitude, but significantly different in phase. The effect of 
the nose absorber on gunner vibration will not be the same as the effect of an 
absorber directly at the gunner station. 


TABLE VII. - VERTICAL 2P (10.8 Hz) VIBRATIONS, g 



Z50 

Nose 

— 

Z90 

Gunner's left 

Z400 

Hori zontal 
stabilizer 

187 kts rolling 
pullout left 

.426/150° 

.464/118° 

.938/68° 

164 kts rolling 
pullout right 

.533/156° 

.494/129° 

.818/61° 

144 kts 

.215/168° 

.201/133° 

.344/44° 

103 kts 45° 
turn 

.121/162° 

.131/116° 

.237/92° 

Sideward flight 
R & L to 35 kts 

.106/152° 

.113/145° 

.249/132° 

Approach and landing 

. 095/ - 69 0 

. 1 00/ -31 ° 

. 271 / -24° 


Figure 1 shows the variation in vertical vibration at the gunner's left 
with variation in the tuning and damping of the nose absorber. The abscissa 
scale is not the same as an RPM sweep because the mobilities and the vibrations 
would change in an RPM sweep. It is impractical to do an RPM sweep in every 
maneuver. However, sensitivity to changes in the tuning frequency are indica- 
tive of bandwidth. 

In most plots of Figure 1, a change in the absorber tuning of less than 1% 
causes a significant change in gunner's seat vibration, indicating an imprac- 
tically narrow bandwidth. Note that 2% structural damping, not zero damping, 
in the absorber gives minimum gunner vibration at 144 knots and in a 103-knot 
turn, but in other maneuvers zero damping gives the minimum. In the rolling 
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Vibration with absorber 
Vibration without absorber 











pullouts 2 % damping reduces the peak of mis tuning more than it raises the depth 
of the valley. The tuning ratio for minimum vibration is very close to 1.0 and 
for gunner's seat attenuation the absorber is tuned to 2P or 10.8 Hz. 

The effect of the nose absorber on the vertical vibration at the hori- 
zontal stabilizer (FS Z400) is shown in Figure 2. The absorber frequency for 
minimum vibration shifts somewhat with maneuvers but is generally lower than 2P. 
Figure 2 also indicates that tuning the nose absorber to 2P (10.8 Hz) to mini- 
mize gunner vibration will increase the horizontal stabilizer vibration by a 
substantial amount in most flight conditions. The converse is true if the nose 
absorber is tuned to minimize horizontal stabilizer vertical vibration. 

There is also a significant change in vibration between zero and 2% damping. 
The extreme sensitivity of vibration at the gunner's station and the horizontal 
stabilizer station to absorber tuning and absorber damping suggests that an 
absorber at the nose would have to weigh much more than 13.61 kg (30 lb) to be 
useful . 

Absorber at the tail rotor . - The principal concern in this example is 
minimization of the lateral vibration at FS Y490, where the vertical stabil- 
izer joins the tail boom. A lateral absorber, weighing 13.61 kg (30 lb), is 
located at the tail rotor gearbox as shown in the sketch on page 30. 
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Vibration with absorber 
Vibration without absorber 



(a) 187 knots rolling 
pullout left. 


.938 

g's 


J 


T 

2 %~> 

\ 

0 % 


(c) 144 knots straight 
and 1 evel . 


.344 

g's 



(b) 164 knots rolling 
pullout right. 



.237 

g's 




Tuned frequency/forci ng frequency 
(e) Sideward flight to (f) Approach and landing. 

35 knots. 


.271 

g's 


Figure 2. - Effect on stabilizer vertical (FS Z400) vibration of 13.61 kg 
(30 lb) vertical absorber at nose (FS Z50) for 0%, 2%, and 
5% absorber structural damping. 







The pertinent acceleration mobilities are 



Y(Y517,Y517) 

< Y(Y490,Y517) - 
Y ( Y440 » Y51 7 ) 


3.597 g/ 1 000 N 
1.576 g/ 1 000 N 
_ .535 g/ 1 000 N 


(1.583 g/100 1 b)/l 7° 
(.693 g/100 1 b )/20° 
(.235 g/100 1 b)/2° 


(27) 


and the flight accelerations are given in Table VIII. 


The lateral response at the absorber attachment point, for 2% absorber 
damping is shown in Figure 3. Comparison of Figure 3 to Figures 1 and 2 (note 
abscissa scale change) shows that the tail rotor gearbox absorber is much less 
sensitive to the frequency ratio, indicating a broader bandwidth than the nose 
absorber. The high lateral driving-point mobility at the tail rotor gearbox 
explains this improvement. As seen from equation (22), the flight accelera- 
tions cancel at the absorber attachment point and Figure 3 is, therefore, the 
same for all flight conditions. 
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TABLE VIII. - LATERAL 2P (10.8 Hz) VIBRATIONS, g 



Y51 7 

Y490 

Y440 

187 kts rolling 
pullout left 

1 .81 5/-91 ° 

1 . 266/ -1 16° 

1 . 282/-1 1 5° 

164 kts rolling 
pullout right 

. 851 /-1 06° 

. 71 0/ -1 21 ° 

. 967/ -1 14° 

144 kts 

1 . 797/-1 33° 

. 878/ —1 50° 

.611/-145 0 

103 kts 45° 
turn 

.405/171° 

.233/175° 

.118/-170 0 

Sideward flight 
R & L to 35 kts 

. 01 7/-45° 

.031 /-54° 

. 07 5/ -43 0 

Approach and landing 

.238/-40 0 

.150/137° 

.184/140° 



Tuned frequency 
Forcing frequency 


Figure 3. - Lateral response at the tail rotor gearbox 
(FS Y517) for 2% structural damping. 
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Comparing Figure 4 with Figures 1 and 2 shows that the tail rotor absorber 
effects on fin lateral vibration are much less sensitive to tuning and absorber 
damping than the nose absorber effects on vertical vibration at the gunner's 
station and stabilizer station. Although the variations in minimum vibration 
with damping are small in Figure 4, the minimum is lowest with zero damping, 
except in the 187-knot rolling pullout where 5% structural damping gives the 
lowest vibration. The absorber tuning frequency for the minimum vibration 
varies somewhat with flight condition but is near 90% of 2P in most cases. In 
this example a tuning frequency of about 9.7 Hz appears to be the best compro- 
mise. Such a selection, as seen from Figure 4, would result in much higher fin 
vibration in approach and landing. Referring to Figure 3, a tuning frequency 
of about 9.7 Hz gives almost the same vibration at the absorber attachment 
point. Figure 5 shows the vibration at FS Y440 for variation in tuning of the 
absorber with 2% absorber structural damping. The 9.7 Hz tuning frequency would 
have a minor effect at this location in most flight conditions. 

The engineer must be cautious of the effects of RPM changes and as a first 
approximation the engineer, assuming the vibrations to remain constant with RPM 
change given no RPM sweep flight data, would utilize the mobility spectrum data 
to create plots similar to those of Figure 4 with an inverted abscissa param- 
eter: i.e., variation in forcing frequency for the tuning frequency selected. 

This would be done on the interactive computer for all locations and directions 
of importance over the RPM range allowable in flight. 
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Vibration with absorber 
Vibration without absorber 


710 

g's 


(a) 187 knots roll ing 
pullout left. 


(b) 164 knots rolling 
pullout right. 


c) 144 knots straight 
and 1 evel . 



mm 


(d) 103 knots 45° turn 


Tuned frequency/forcing frequency 
(e) Sideward flight to (f) Approach and landing. 

35 knots 


Figure 4. - Effect on fin lateral (FS Y490) of 13.61 kg (30 lb) lateral 

absorber at tail rotor for 0%, 2 %, and 5% absorber structural 
damping. 









Vibration with absorber 
Vibration without absorber 


Figure 5. 



1.282 
g 1 s 


(a) 187 knots rolling 
pullout left. 



and level . 


2i 


C 



.5 


1.0 


1.5 



.967 

g's 


(b) 164 knots rolling 
pullout right. 



.118 

g's 


(d) 103 knots 45° turn 



.184 

g's 


1.0 


1.5 


(e) Sideward flight to 
35 knots. 


Tuned frequency/forcing frequency 

(f) Approach and landing 


- Effect on boom lateral (FS Y440) of 13.61 kg (30 lb) lateral 
absorber at tail rotor for 2 % absorber structural damping. 
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Active Vibration Suppression 


Vibration suppression and vibration isolation are distinguished by con- 
sidering suppression to imply something other than separating (isolating) the 
externally excited structure. Active vibration suppression devices may be 
applied without consideration of the physical locations of the external exciting 
forces and moments. In this report the term active is used in the customary 
sense to indicate a powered device which creates an external force with magni- 
tude and phase controlled by reference to a feedback signal of vibration or 
strain. The vibration or strain to be controlled will generally not be at the 
location of the active vibration suppressor. 


Active vibration suppression with time-domain control has been highly 
effective in some applications at frequencies far below the lowest blade passage 
frequency in helicopters and continuing work in this specialized area may even- 
tually be important at helicopter frequencies, but such matters are beyond the 
scope of this report. Frequency-domain control of active vibration suppressors, 
perhaps by minicomputers, is not as frequency limited as time-domain control 
and is close to the state-of-the-art of shaking hardware. For these reasons, 
frequency-domain control is implied in the considerations in this report. It 
is to be understood in the discussions following that any active vibration 
suppressor can, by the methods described for a given frequency, be simultan- 
eously applied to other frequencies which are not necessarily limited to har- 
monics of main rotor blade passage frequency. 


Minimizing the norm of a set of flight accelerations or strains . - Active 
vibration suppressors, R in number, can be used to minimize the sum of the 
squares of K, greater than R, accelerations or strains. Let the f complex 
vector be that of the forces or moments applied by the R active vibration sup- 
pressors and the primed quantities be the resulting complex accelerations and 
strains. The matrix equation which defines this situation is 


f 




' „ 

3 . 


^a 


Y ar 


. — - 

— 

h + 

— 

e s 


E s 







sr 


K x 1 K x 1 K x R 



R x 1 


( 28 ) 
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or, more concisely. 




kr 


K x 1 K x 1 KxR R x 1 


(29) 


The sum of the squares (Euclidian norm) becomes 




f ’ 

l*T / 

Kj 

f V 

M 

i*t r 

l rj 

_ 


kr 


*T 


*T 


kr 


kr 


f r + f r 


*T 


kr 


*T 


(30) 


where * denotes the complex conjugate. 


The left hand side of equation (30) is, of course, a real scalar and the 
minimum sum of the squares is obtained by setting the partial derivative of this 
complex scalar function with respect to the complex transpose of the f vector 
to zero. 





(31) 


The complex forces from the R active vibration suppressors necessary to mini- 
mize the sum of the squares (the Euclidian norm or sum of the squares of the 
absolute values) of the K (greater than R) vibrations or strains in any flight 
condition are 


r \ 



This process is, except for the minus sign, identical to that of Force Deter- 

. .. 5 

mmation. 
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If R were equal to K, then all K accelerations and strains would be zero. 
In this case it is necessary to examine the vibration or strain of each impor- 
tant motion coordinate which is not of the set of K = R because some of those 
of the nonnulled set could have been made much worse. When K is greater than 
R, it is necessary to check every significant motion coordinate, including 
those of the minimized set, because a reduction in the sum of the squares of 
the minimized set does not guarantee that some motion coordinates in the mini- 
mized set are not amplified beyond tolerable values. The changed vibration is 
determined from 


q '• = q . + 
lxl 1x1 



( 34 ) 


Some Types of Active Vibration Suppressors 


The following types of active vibration suppressors are discussed to illus- 
trate the application of analytical testing to active vibration suppression. 

The list is not intended to be all inclusive and only a cursory examination of 
the practicality of the devices is presented. 


Active mass . - The shaking of a reaction mass by an actuator attached to 
the airframe creates an external force on the airframe as shown below. The 
response, acceleration or strain, to be controlled will usually be remote from 
the active mass and may be in any spatial direction. The feedback to the 


ACCELEROMETER 
OR STRAIN GAGE 



POWER 

SOURCE 
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controller from the shaker may be both force magnitude and phase. To obtain 
the desired force level on the airframe with minimum reaction mass weight 
requires the reaction mass to move through a large displacement. 

Active stiffness . - A shaker or actuator may be placed between two points 
on the structure as shown in the following sketch. 



The response at j is given by 


q j q j + (V jx2 ' Y jxl } f x 


(35) 


where xl and x2 must be col inear. Feedback from the shaker may be both force 
magnitude and force phase. 


Control surface suppressor . - An aerodynamic surface, such as a canard or 
horizontal stabilizer, can be vibrated as shown in the following sketch to 
create an external force. Considering the control surface to be essentially 
rigid, one may obtain phase feedback. However, the actual external vibratory 
force produced may not be known accurately unless the suppressor can be test 
calibrated with a suitably placed strain gage or bending bridge because of the 
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difficulties in calculation of the force from unsteady aerodynamics and aero- 
dynamic interferences. 



A disadvantage of using a fixed-system control surface as a vibration sup- 
pressor is airspeed dependence because the maximum force required may not be at 
the highest airspeed. 

Rotor control suppression . - Variation in rotor pitch is another possible 
method of active vibration control. A schematic for vibratory tail rotor pitch 
is shown below. 

ACCELEROMETER 



In some types of tail rotor installations the vibratory force magnitude 
and phase may be obtained from strain gages on the tail rotor shaft. 
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Active absorber . - The active absorber can be used to suppress vibration 
at a particular frequency. The major advantage of an active absorber is the 
reduction of the shaker size required for a given output force. Figure 6 
depicts a conventional absorber with damping. The active absorber is simply 
a spring in parallel with the shaker, as shown in Figure 7. 


MASS 



Figure 6. - Conventional 
absorber with damping. 



Figure 7. - Active absorber. 
Damping may be zero. 


The displacements of the absorber mass and the fuselage are obtained from 


q 


a 


q 


n 



Y' 


ra 



(36) 


where the elements of the matrix in equation (36) are displacement mobilities 
with the absorber on the fuselage. Substituting for the absorber parameters 
leads to the result 


where 


' Y' 
aa ar 

_ 1 

K(l+ig) + Y^ 1 

K( 1+ig) 

' Y' 
ra rr 

A 

K(l+ig) 

K(l+ig) - oo 2 m 


A = K(1+ig) (Y r ~ 1 - Jm) - Jm V^ 1 


(37) 


(38) 


and Y^ r is the driving-point mobility without the absorber. The spring stiff- 
ness, spring structural damping, and absorber mass are K, g, and m, respectively. 
The motion of the absorber mass is obtained from equation (36) and, after elim- 
inating f , becomes 
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( 39 ) 



For a shaker force f on the absorber mass and -f on the fuselage, equation (28) 
determines the displacement of the absorber mass as 


q 


a 





(40) 


and the displacement of the fuselage as 


q 


r 





f 


then 

fy 1 

j 



q 1 1 - q ' 1 = 
M a H r 

ar - 1 

q 1 + 

Y 1 + Y' 

- 2 Y 1 

Vr J 


aa rr arj 


f 


(41) 


(42) 


Substituting for the driving-point displacement mobility of the absorber, equa- 
tion (22) leads to the result 


q r q r 


[ K(l + iq) - A] 

[K(l+1g) - w 2 m - Y rr K(l + ig) u 2 m] 


(43) 


where q^ is the motion of the fuselage without the absorber. 


Equation (42) can then be written as 


q 


a 


w 2 mq r + (1 - w 2 m Y rr ) f 
K(l+ig) (1 - w 2 m Y rr ) - w 2 m 


(44) 


The force on the fuselage is 

f r = K(l+ig) (qj' - q; 1 ) - f 


(45) 
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or 


2 

oj m 


f + q r K(l + ig )] 


r K(1 + i g ) ( 1 - Y rr w 2 m) - i/m 


(46) 


and in terms of acceleration 


f = 

r 


- iynd + ig) + ^2 f 


1 - 


+ m 


f Y R - gY 11 


+ im 


fy I + gY ^ 

rr y rr 


+ ig 


(47) 


where S2y is the antiresonant frequency at r created by the absorber, 


To minimize the denominator with tuning, let 


= m 


Y R - gY ^ 
rr a rr 


+ 1 


(48) 


which sets the real part of the denominator of equation (47) to zero. The 
tuning of the active absorber is independent of flight vibrations at the attach- 
ment point and the response point whereas the tuning of a conventional absorber 
is not. With the tuning condition of equation (48), equation (47) becomes 


f = i 
r 


q r m (1 + ig) - f w 2 /fi T 2 


^Y 1 + gY R 
rr y rr 


(49) 


+ g 


From equation (48) the tuning frequency of the active absorber is 


fiy = 


(50) 


/ 


+ m 


Y R - gY ^ 
rr y rr 


In operation, the minicomputer controller determines the magnitude and phase 
required of the shaker from 

i 


f = 


fin 


VJ + 9 ? rJ] + Wf + 90< 


+ lv m 


2 ) -1 
1 + 9 / 4> r + tan g 


(51) 
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The force on the fuselage at r, f , for zero vibration at j is determined from 


' Y. 

jr 

There are two feedback signals to the controller of the active absorber: the 

acceleration at the attachment point and the vibration (strain or acceleration) 
of the motion coordinate to be nulled. 


Examples of Active Vibration Suppression 

The applicability of analytical testing for examining the effects of 
active vibration suppressors on airframe vibration is illustrated using AH-1G 
ground and flight test vibration data. Equations (33) and (34) are used to 
determine the required control forces and the changed vibrations. 

Horizontal stabilizer aerodynamic suppressor . - This example considers the 
effects on vertical vibration for 2P excitation of the horizontal stabilizer. 

As shown in the following sketch, the horizontal stabilizer is vibrated in pitch 
with sufficient force magnitude and phase to give zero vertical vibration at the 
pilot's seat (FS Z140). The effects at FS Z90, FS Z400, and FS Z540 are also 
examined. The vertical acceleration mobilities at these four locations are 
shown in Table IX. Table X presents the flight vertical accelerations for four 
flight conditions. 
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Z540 



At 187 knots, the required force is 

f = . q(Z140) = _ . 1 47q 1 s/1 24° 

r Y(Z400,Z140) ( .065g 1 s/1000 N)/-84° 


(53) 


- - 2261 N ( -51 0 1 b)/208° 


and similarly for other flight conditions. The vibration at any point is 
obtained from equation (34) as 


q '. = q . + Y • f 
Jr r 


(54) 


and the results are shown in Table XI. 
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TABLE IX. - ACCELERATION MOBILITIES AT 10.8 HZ, g/1000 N (g/100 lb) 



Z90 

Z140 

Z400 

Z540 

Z90 

.103/10° 

(.046) 

.070/6° 

(.031) 

,038/-138° 

(.017) 

.288/9° 

(.128) 

Z140 

.070/6° 

(.031) 

.052/5° 

(.023) 

O 

oo 

un cri 

CO CM 

o o 

.124/25° 

(.055) 

Z400 

•038/-138 0 

(.017) 

.065/-84 0 

(.029) 

.072/64° 

(.320) 

. 672/-32 0 
(.299) 

Z540 

.288/9° 

(-128) 

.124/25° 

(.055) 

. 672/- 32 ° 
(.299) 

2.855/-8 0 

(1.270) 


TABLE X. - FLIGHT ACCELERATIONS AT 10.8 HZ, g 



187 Knots 
rol 1 ing 
pul 1 out 
left 

164 Knots 
rol 1 ing 
pullout 
right 

144 Knots 
straight 
and level 

103 Knots 
45° turn 

Z90 

.335/120° 

.322/134° 

.118/120° 

.064/124° 

Z140 

.147/124° 

.274/114° 

.114/99° 

.078/95° 

Z400 

.938/68° 

.818/61° 

.344/44° 

.237/92° 

Z540 

1 . 992/-1 18° 

1 . 454/ - 1 31 ° 

. 769/ -1 31 ° 

.684/-1 67° 


Assume that the horizontal stabilizer is vibrated in pitch with sufficient force 
magnitude and phase to give zero vertical vibration at the pilot's seat, as 
shown by the previous sketch. 
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TABLE XI. VERTICAL VIBRATION AT 10.8 HZ WITH AND WITHOUT 
. HORIZONTAL STABILIZER AERODYNAMIC SUPPRESSOR, g 



187 Knots 
roll ing 
pullout 
left 

164 Knots 
rolling 
pul lout 
right 



144 Knots 
straight 
and level 

103 Knots 
45° turn 

Force at 
Z400 

-2255 N/208° 
(-507 lb) 

-4204 N/l 98° 
(-945 lb) 

-1748 N/l 83° 
(-393 lb) 

-1197 N/l 79° 
(-269 lb) 


With Without 

With Without 

With Without 

With Without 

Z 90 
Z140 
Z400 
Z540 

.287 .335 
0 .147 
2.51 .938 
1.95 1.992 

.316 .322 
0 .274 
3.798 .818 
2.523 1.454 

.120 .118 

0 .114 

1.58 .344 

1.263 .769 

.074 .064 

0 .078 

1.074 .237 

.592 .684 


The pilot's seat vibration is zero in all flight conditions with a large 
increase in vertical vibration at the horizontal stabilizer station on the boom. 
The forces required are very large, however. 

For purposes of illustration let it be assumed that the horizontal sta- 
bilizer has an area of approximately 1 square meter and that trim requirements 
would permit a 2P pitch vibration of + 3° maximum producing a vertical force of 
1112 N (250 lb) at 187 knots airspeed. With the maximum force proportional to 
the square of the airspeed, the vibrations obtainable with this arrangement are 
given in Table XII. 

Even with the very large forces of Table XI, zero vibration at the pilot's 
seat from an active vibration suppressor at the horizontal stabilizer station 
is obtained at the expense of large increases in tail boom vibration with neg- 
ligible changes in gunner's seat vibration. The reduction in vibration at 
the pilot's seat and gunner's seat using the horizontal stabilizer as an aero- 
dynamic suppressor, as shown in Table XII, are not impressive. 
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TABLE XII. - HORIZONTAL STABILIZER FORCED AT 10.8 HZ TO 
MINIMIZE PILOT'S SEAT VIBRATION 



187 knots 
rol 1 ing 
pullout 
left 

164 knots 
rol 1 i ng 
pullout 
right 

144 knots 
straight 
and level 

103 

45° 

knots 

turn 

Force at 
FS Z400 

-1112 N/208° 
(-250 lb) 

-854 N/l 98° 
(-192 lb) 

-658 N/183° 
(-148 lb) 

-338 N/l 79° 
(-76 lb) 

Vertical vibration with and 
without suppression, g 


Wi th 

Without 

With Without 

With 

Without 

With 

Without 

Z90 

.309 

.335 

.315 .322 

.114 

.118 

.064 

.064 

Z140 

.075 

.147 

.218 .274 

.071 

.114 

.056 

.078 

Z400 

1.700 

.938 

1.409 .818 

.802 

.344 

.465 

.237 

Z540 

1.821 

1.992 

1.298 1.454 

.804 

.769 

.551 

.684 


T-tail aerodynamic suppressor . - The effects on vertical vibration for 2P 
excitation of the T-tail, as shown below, are illustrated in this example. 



The principal objective is to give zero vertical vibration at the gunner's seat 
(FS Z90) . The airframe locations and flight conditions are identical to the 
previous example. The T-tail horizontal control surface is not required for 
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trim and can be operated at higher 2P vibratory angles of incidence than the 
horizontal stabilizer (FS Z400). Therefore, assuming a T-tail area of approxi- 
mately 0.4 square meter, a vibratory force of about 1200 N (270 lb) at 164 knots 
can be generated. 

The results of using the T-tail aerodynamic suppressor to null the gunners' 
seat vibration are shown in Table XIII. With the 164 knot presumption, the 
required forces are well within the airspeed-squared requirement. The T-tail 
is outside the path of the main rotor tip vortex indicating that it may not be 
the source of increased external loads at main rotor 2P. 

Significant vibration reductions result at the pilot's seat and tail boom 
with zero vibration at the gunner's seat. Whether the increase in vibration 
at the tip of the fin (FS Z540) is structurally significant and whether the 
lateral offset of the T-tail force is detrimental or beneficial to fin bending 
are questions that would be answered by further analsis using analytical 
testing. 


TABLE XIII. - VERTICAL VIBRATION AT 10.8 HZ WITH AND 
WITHOUT T-TAIL AERODYNAMIC SUPPRESSOR, g 



— 

187 Knots 
rolling 
pullout 
left 

164 Knots 
rolling 
pullout 
right 

144 Knots 
straight 
and level 

103 

45° 

Knots 

turn 

Force at 
FS Z540 

-1164 N/l 11° 
(-261.7 lb) 

-1119 N/l 25° 
(-251.6 lb) 

-410 

(-92. 

N/l 11° 
2 lb) 

-222 

(- 

.4 N/l 1 5° 
50 lb) 


With Without 

With Without 

Wi th 

Wi thout 

With 

Without 

Z 90 

0 .335 

0 .322 


.118 

0 

.064 

Z140 

.031 .147 

.181 .274 

.080 

.114 

.062 

.078 

Z400 

.226 .938 

.437 .818 

.197 

.344 

.092 

.237 

Z540 



5.000 1.992 

3.975 1.454 

1.740 

.769 

1.295 

.684 
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Active absorber at the gunner's seat . - An active vertical absorber at the 
gunner's seat (FS Z90) is controlled to provide zero vertical vibration at the 
pilot's seat (FS Z140). The effects on vibrations were calculated at the loca- 
tions shown in the sketch below for four flight conditions. The required 
acceleration mobilities and flight accelerations are given in Tables IX and X, 
respectively. 



The results presented in Table XIV indicate that the pilot's seat (FS Z140) 
vertical vibration is zero for all flight conditions. The gunner's seat 
(FS Z90) vertical vibration remains the same in the 103-knot turn but is 
reduced in the other maneuvers. The horizontal stabilizer (FS Z400) vertical 
vibration does not change appreciably and the fin (FS Z540) vertical vibration 
is increased. 

The active absorber parameters are calculated from equations (50), (51), 
and (52). At 187 knots the suppressor force is -2100 N (-472 lb)/118° acting 
through .335 g/120° at 10.8 Hz. The output power is the real part of the pro- 
duct of force and velocity which gives 33.08 Nm/s (.044 hp) consumed in the 
suppression. For a 4.536 kg (10 lb) reaction mass, the required tuning fre- 
quency is 10.78 Hz and the required shaker force is 65.5 N (14.8 lb)/41.1°, 
assuming 3 percent hystegetic damping in the absorber. 


49 


TABLE XIV. - VERTICAL VIBRATION AT 10.8 HZ WITH AND WITHOUT 
ACTIVE ABSORBER AT GUNNER'S STATION, g 



187 Knots 
rolling 
pullout 
left 

164 Knots 
rolling 
pullout 
right 

144 Knots 
straight 
and level 

103 

45° 

Knots 

turn 

Force at 
FS Z90 

-2108 N/118° 
(-474 lb) 

-3732 N/l 08° 
(-839 lb) 

-1637 N/93° 
(-368 lb) 

-1121 N/89° 
(-252 lb) 


With Without 

With Without 

With 

Without 

With 

Without 

Z90 

.123 .335 

.132 .322 

.066 

.118 

.064 

.064 

Z140 

.000 .147 

.000 .274 

.000 

.114 

.000 

.078 

Z400 

.939 .938 

.834 .818 

.349 

.344 

.272 

.237 

Z540 

2.315 1.992 

2.151 1.454 

1.117 

.769 

.781 

.684 









Active absorber at the tail rotor . - In this example the objective is to 
produce zero lateral vibration on the fin (FS Y490) at 10.8 Hz. An active 
lateral absorber at the tail rotor gearbox was selected as the vibration sup- 
pressor, although excitation of tail rotor pitch is an alternative. The pre- 
ference of harmonic control of the tail rotor or an active absorber depends on 
many factors among which are tail rotor blade loads. Main rotor 2P excitation 
of the tail rotor may increase or decrease tail rotor fatigue life and this 
requires additional investigation. The effects on lateral flight accelerations 
were calculated at the locations shown in the following sketch for four flight 
conditions. The required acceleration mobilities and flight accelerations are 
shown in Tables XV and XVI, respectively. The flight conditions are identical 
to the previous example. 

The results shown in Table XVII indicate that the tail rotor gearbox vibra- 
tion is negligibly affected in the rolling pullouts, but significantly reduced 
in 144- knot level flight and in the 103-knot turn. The vibration at FS Y440 is 
reduced substantially in the 187-knot left rolling pullout and at 144-knot level 
flight, reduced somewhat in the 164-knot right rolling pullout, and increased in 
the 103-knot turn. Changes in pilot and gunner vibration are not significant. 

As expected, the fin (FS Y490) vibration is zero for all flight conditions. A 
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TABLE XV. - ACCELERATION MOBILITIES AT 10.8 HZ RELATIVE TO 

FS Y51 7 (TAIL ROTOR LATERAL), g/1000 N (g/100 lb) 


FuseTage stations 


Y51 7 

Y490 

Y440 

Y140 

Y90 

3.597/17° 

1.576/20° 

.535/2° 

.178/103° 

.339/84° 

(1.600) 

(.701) 

(.238) 

(.079) 

(.151) 


TABLE XVI. - LATERAL FLIGHT VIBRATION, g 


Airspeed 

Fuselage stations 

Y 51 7 

Y490 

Y440 

Y140 

Y90 

187 kts 

1 . 81 5/-91 ° 

1 . 266/— 1 1 6° 

1 . 282/-1 1 5° 

. 1 23/ -96 0 

. 1 91 /-94° 

164 kts 

.851/-106 0 

. 71 0/-1 21 ° 

. 967/ - 1 1 4° 

. 1 1 8/-89° 

. 1 47/-83° 

144 kts 

1 . 637/-1 38° 

. 87 8/ - 1 50° 

.61 1 /-I 45° 

.055/-1 14° 

. 1 03/-1 20° 

103 kts 

.405/171° 

.233/175° 

. 048/-95 0 

. 048/-95 0 

.072/-105 0 
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TABLE XVII. - LATERAL VIBRATION AT 10.8 HZ WITH AND WITHOUT 
ACTIVE ABSORBER AT TAIL ROTOR GEARBOX, g 



187 Knots 
rol ling 
pullout 
left 

164 Knots 
rol 1 ing 
pullout 
right 

144 Knots 
straight 
and level 

103 

45° 

Knots 

turn 

Force at 
FS Y 51 7 

-803 N/-1 36° 
(-181 lb) 

-450 N/-141 ° 
(-101 lb) 

-557 N/-1 70° 
(-125 lb) 

-148 N/l 55° 
(-33 lb) 


With Without 

With Without 

With 

Without 

With 

Without 

Y51 7 

1.552 1.815 

.857 .851 

.602 

1.657 

.127 

.405 

Y490 

0 1.266 

0 .710 

0 

.878 

0 

.233 

Y440 

.887 1.282 

.755 .967 

.356 

.611 

.104 

.048 

Y140 

.140 .123 

.092 .118 

.073 

.055 

.022 

.048 

Y 90 

.183 .191 

.068 .147 

.119 

.103 

.027 

.072 


comparison of the tail boom (FS Y440) vibration reduction given in Table XVII 
to Figure 5 shows that the active absorber behaves quite differently from the 
conventional absorber at stations other than that suppressed. 

The active absorber parameters are calculated from equations (50), (51), 
and (52). At 187 knots the suppressor force is -803 N (-181 1 b )/-l 36° acting 
through 1.815 g/-91° at 10.8 Hz which gives 172 Nm/s (.229 hp) consumed in the 
suppression. For a 2.27 kg (5.0 lb) reaction mass, the required tuning fre- 
quency is 10.41 Hz and the required shaker force is 38.9 N (8.76 lb)/174°, 
assuming 5 percent hysteretic damping in the absorber. 


Stiffness Changes 


Skin and plate changes . - The basic 
obtained by combining equations (11) and 


equation for a stiffness change is 
( 12 ) 


-2L 


jk_ 


V 

, 

k. — * 


— 

yl'cT 

. £ j . 


, . 


Y kr 

L_ 4 




( 55 ) 
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where p< rr j is the stiffness matrix of the skin section as a free body, {q } is 
the vector of flight displacements at the change coordinates and |~Y 1 is the 
matrix of displacement mobilities of the airframe at the change coordinates. A 
skin, as opposed to a plate, has negligible transverse stiffness and the coor- 
dinates, |q r |, are on the surface of the airframe. The differences among the 
coordinates, jqL are most important in a stiffness change. It would be 
impractical in most cases to obtain jq r j from flight accelerometers alone 
because of the small difference between large numbers. This problem, however, 
is resolved using strain gages. 

Figure 8 shows a skin section and the rectangular coordinates of the rele- 
vant stiffness matrix. The skin may be of nonhomogeneous materials with noncon- 
stant thickness, such as a two dimensionally tapered composite, with directional 
modulii of elasticity or a simple metal sheet. Only the stiffness matrix needs 
to be known in applying equation (55). The effects of proportional or nonpro- 
portional damping may be included in equation (55) by replacing |~K 1 with a 
complex matrix which is the sum of the real stiffness matrix and an imaginary 
damping matrix, frequency dependent for viscous damping and frequency indepen- 
dent for hysteretic damping. 

It is immaterial to the analytical testing approach whether the surface of 
the helicopter at which the skin or plate change is to be made is flat or 
curved. 

For n nodes on a skin there are 2n rectangular coordinates and the stiff- 
ness matrix (the matrix of the partial derivatives of force to displacement) is 
a square matrix of order 2n and rank 2n - 3. For practical flight and shake 
test measurements the skin section may be modeled in terms of 2n - 3 strain 
coordinates as shown in Figure 9. 

The displacement mobility matrix for the strain coordinates of, for example, 
the skin section of Figure 9 may be determined as shown in Figure 10 where equal 
and opposite forces are applied across a distance 6. This type of displacement 
mobility (influence coefficient) matrix becomes 


53 




Y (C) . 

L J 3Af . 

J 


(56) 


which is, of course, nonsingular and may be inverted to give a stiffness matrix 

, 1-1 


■(c) 


= [k] = 


3Af - 


9Aq, 


(57) 


where 


and 


Af j ' i (f j + i - f j-i > 


Aq k q k+l " q k-l 


If, on the other hand, the conventional type 2n x 2n stiffness matrix K 
is available the strain type stiffness matrix k may be found using a coordinate 
transformation. Note that 


[K] = 



(58) 
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A transformation matrix [T] exists such that 


where 


Let 


{Aq} = [T] {q} 
(2n-3}xl (2n-3)x2n 2nx1 


[T] 


W 

3q 


2n 



1 


3 (2n-3) 


( 59 ) 


(60) 


be the eigenvectors of K where the superscript R indicates the rigid body modes, 
those of zero eigenvalues. 

[K] [*] = w [A] (61) 

and 

M T [K] [$] = CK *3 


for nonzero eigenvalues. The inverse of K may be expressed as 


[K] 


-1 



7=r ! 0 


r -\" 

0 jk 



(62) 


This matrix is singular because of the three reciprocal of zero terms in the 
diagonal matrix. The 2n x 2n eigenvector matrix is necessarily nonsingular. A 
coordinate transformation of equation (62) using equation (60) yields 


[k] 


-1 


9 Aq 
dhf 


[T] 


o 

0 Ik 



(63) 


which is a (2n-3) x (2n-3) matrix. The transformation of equation (60) applied 
to a rigid body mode is necessarily zero. Therefore, 
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[kj 


0 ! T$ 


L i 


- n 5 :» 

->T 
0 1 — 
u 1 K* 


0 ' T* 


(64) 


The 0 x 0/0 indeterminant terms in the product are zero and equation (64) may 
be written as 


r(c) 

rr 


= [k]’ 1 = [T]M 




[*] T [T] T 


from which it follows that k, a nonsingular matrix of order (2n-3), is given by 


3 Af" 


LBAqJ 


= [k] = | [T][4.] r T fK*] [T]0] 
= |[T][t]|' T [*] t [k][>] 


-1 


[T][*] 


1-1 


(65) 


It is also necessary to establish the types of mobility measurements which 
are practical in an actual airframe. 


For forces f. + ^ and f j _ -j only in. Figure 9 


4 + fVi , ^k-i . 

Aq k - »f~ Vl aq.f f j-l - 3f j+1 - f j + l 


3 Vl , 
3 f j-l 


( 66 ) 


But, -f. , = f.,, = Af.. Therefore 
J-l J+i J 


^k 

Af j 


3q 


k+1 3q k+l 


3q 


k -1 + 3 q k-l 


3 Vi 9 f j-i 3 f j+i 3 f j-i 


N 

= I 

i=l 


u/ \u _ u; u/ -U/ vi/ 

Y k+l,i Y j+l,i K+l,i ^ j -1 , i k -1 ,i T j +1 ,i 


+ \-l ,i Y j -1 , i I F i 


N 

" ih 


\+l ,i ~ ^k-l ,i 


¥ j+l J " >i 


F i 


(67) 
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For forces only at j 



f 


' -Y 

Aq. N 

* = y 

Vl,i " Vl,i 


Vi.i - *1-1.1 

Af. .A 

J 1 = 1 

S k J 


*7 1 


6 j \ F i 


( 68 ) 


or, for any number of forces. 


[ Y rr] 




3Af 


7, * ( ki ’ ( ji s k F i 


(69) 


Simi larly. 


3q. N / x 

iAf; ■ I "ki 5 r F i 


(70) 


With a strain stiffness matrix, equation (55) becomes 

il-l 






- 


k 

_ 

ik 


3 V Af r 


i 

p „ 




3e i 




, e j. 


3Af r 



ru + r k i 

* l rr 


3Aq, 


3 Af 


P rr 


£ r 6 r 


(71) 


In the case of a plate, there is transverse stiffness and nonnegl igible 
mass. The aircraft is tested with n accelerometers perpendicular to the sur- 
face, for n attachment points, in addition to the strain gages, as shown in 
Figure 11 . 

Numerical example of a transformation to a strain stiffness . - The six 
pinned rods of Figure 12 provide a simple illustration of the transformation 
from a stiffness matrix with rectangular coordinates to one with strain coordi- 
nates . 

In the stiffness matrix with rectangular coordinates, shown following 
Figures 11 and 12, the k terms are axial spring rates subscripted with the ter- 
minal positions of the rods. 
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[K] 


2 

k 12 +k 13 cos “ 

-k 12 

-ki3 c °s 2 ° 

0 

" k 12 

k 2 ^cos a+k^ 2 

0 

-k 24 cos Z a 

-k 13 cos 2 a 

0 

k 34 +k 13 cos a 

~ k 34 

0 

2 

-k^CGS a 

* k 34 

k 24 c °S 2 a+ k 34 

k^SfnaCOSa 

0 

-k^sJnaCOSa 

0 

0 

-k^SlnaCOSa 

0 

k 24 SlnaCOSa 

-k^ ^sinaCOSa 

0 

k| 3 sinaCOSa 

0 

0 

k^slna 

0 

-k 24 SlnaCOSa 

X 1 

x 2 

x 3 

x 4 


k 13 cosasino 0 

-k^SinaCOSa 0 

0 -k 24 COSaSlna 

0 k 24 $1naCOSa 

■k^sinacosa 0 

k^SinaCOSa 0 

0 k 24 SinaCOSa 

0 -k 24 SinaCOSa 

2 

k 14 +k 13 s * n a 0 

-k 13 sin 2 a -k 14 

2 

0 k 23 +k 24 s ’ n a 

. 2 

"*^23 —k^^sin & 

' k 13 sin a ” k 23 

2 

^23 + * c 13 S ^ n ^ ^ 

. 2 

“ k 14 *k24 sl(1 a 

0 k 14 +k 24 sin2a ' 

2 ) z 2 

z 3 z 4 


( 72 ) 


This matrix has a degeneracy of three since the null vector is equal to the 
product of the K matrix and the transpose of the following independent vectors: 


1 


representing x translation, 

|0 0 0 0 1 1 1 1 

representing z translation and 

|0 0 -sina -sina 0 cosa cosa 0 

representing a rotation. 
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The coordinate transformation equation in this case is 














x i 

X 1 " 

X 2 


1 

-1 

0 

0 

0 

0 

0 

0 


x 2 

X 4 " 

x 3 


0 

0 

-1 

1 

0 

0 

0 

0 


X 
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x 4 

' Z 1 - 

Z 4 


0 

0 

0 

0 

1 

0 

0 

-1 

. 

> 

Z 1 

Z 2 

z 3 
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0 

0 

0 

0 

1 

-1 

0 


Z 2 

D 1 - 

°3 


COSa 

0 

-COSa 

0 

sina 

0 

-si na 

0 


Z 3 

. 

> 











Z 4 
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or 


T Aq } = [T] {q} 

Let the spring rates, force/deflection, in Figure 12 be 


[K] = 



k l 2 

k 34 = 

3 x TO 5 






k l 4 

k 23 = 

2 x 10 5 






k l 3 

k 24 = 

5 x 10 5 





i ffness 

matrix, i< 

, i n the 

X and Z 

coordinates is 



700000 

-300000 

-400000 

0 

200000 

0 

-200000 

0 

-300000 

700000 

0 

-400000 

0 

-200000 

0 

200000 

-400000 

0 

700000 

-300000 

-200000 

0 

200000 

0 

0 

-400000 

-300000 

700000 

0 

200000 

0 

-200000 

200000 

0 

-200000 

0 

300000 

0 

-100000 

-200000 

0 

-200000 

0 

200000 

0 

300000 

-200000 

-100000 

-200000 

0 

200000 

0 

-100000 

-200000 

300000 

0 

0 

200000 

0 

-200000 

-200000 

-100000 

0 

300000 
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The transformation matrix is 



"1 

-1 

0 

0 

0 

0 

0 

O' 


0 

0 

-1 

1 

0 

0 

0 

0 

[T] - 

0 

0 

0 

0 

1 

0 

0 

-1 


0 

0 

0 

0 

0 

1 

-1 

0 


894427 

0 

-.894427 

0 

.447214 

0 

-.447214 

0 


The eigenvectors of nonzero eigenvalues of K are 



"-1 

-1 

1 

/2-1 

O' 


1 

-1 

-1 

1-/2 

0 


1 

1 

1 

1-/2 

0 


-1 

1 

-1 

/2-1 

0 

M = 

1-/2 

-1/2 

0 

-1 

-1 


1-/2 

1/2 

0 

-1 

1 


/2-1 

1/2 

0 

1 

-1 


_/2-l 

-1/2 

0 

1 

1 
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(76) 


The eigenvectors of K are orthogonal because K is symmetrical. The strain 

stiffness matrix, k, which is the matrix of 3Aq./3Af., is found using equa- 

K J 

tion (65). 


675781.2 

375781.2 

187890.6 

187890.6 

-419262.7 

375781.2 

675781 .3 

187890.6 

187890.6 

-419262.7 

187890.6 

187890.6 

293945.3 

293945.3 

-209631 .4 

187890.6 

187890.6 

293945.3 

293945.3 

-209631 .4 

-419262.7 

-419262.7 

-209631 .4 

-209631 .4 

1E+06 


The inverse of the above k matrix is a strain mobility matrix which is the 
same as a strain influence coefficient matrix and is shown on the following 
page. 
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2.54902E-06 


-7.843137E-07 


-5.882353E-07 


4.932503E-07 


Ck] 


-1 


-5.882353E-07 


-7.8431 37E-07 

2.54902E-06 

-5.882353E-07 

-5.882353E-07 

4.932503E-07 

-5.882353E-07 

-5.882353E-07 

4.558823E-06 

-4.411 764E-07 

3. 699378E-07 

-5. 882353E-07 

-5.88Z353E-07 

-4.411764E-07 

4.558823E-06 

3.699377E-07 

4.932503E-07 

4.932503E-07 

3.699378E -07 

3.699377E-07 

1.568704E-06 


(78) 


As shown by this simple illustration, the coordinates conventionally used 
in finite element analyses for skin stiffness changes of any number of nodes 
can be changed by matrix transformations to accommodate strain instrumentation 
only, a matter of critical importance in practice where differential col inear 
displacements along a fuselage surface are very small. 


Application of analytical testing to a skin change . - From equation (71), 
it is seen that the flight accelerations with the skin change are 



3q k /3A f r 



r i 


~3Aq ' 

m + 

V. 


r 

3Af r 



(79) 


and the flight strains with the change are 







3Aq j/ 3Af ^ 


M + 

l<rr 






~3Aq r 


3Af 


L 




Defi ne 



m + 



_ 3Aq r 


-1 

k rr 

3Af _ 



r 



L_ J 







6 


(80) 


(81) 


For analytical purposes, this may be thought of as a new set of external forces 
operating on the aircraft and caused by the skin change. In this case, these 
are tensile forces, as illustrated in Figure 10. Equations (79) and (80) may be 
written 
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I 



The strain stiffness matrix, k rr , of the contemplated skin change is obtained 
from a finite element analysis as discussed above. From flight tests of the 
baseline aircraft, the accelerations {q . } and the strains { e . } at coordinates 
of interest are known. From modal shake testing (or finite element analysis) 

of the baseline aircraft the orthonormal modes [y,] and and the fre- 

quency function [F] are known. 

The nonflying shake test aircraft is instrumented with strain gages as 
shown in Figure 13 for a nine node change and a modal analysis test is done to 

r fell 

obtain the orthonormal modal elements of ' . The distance between nodes 
should be as large as possible to avoid large orthonormal strain mode elements 
at very high frequencies which occur when there is negligible inertial effect 
at a node. Long gage strain gages are preferred to minimize local effects. 



Figure 13. - Strain gages on the fuselage for flight and shake tests. 


64 



In equation (82) 



m 



and in equation (83) 





(84) 


(85) 


A certain amount of engineering judgment must be applied to the mobilities 
and stiffnesses used in the equations when dealing with high load factor maneu- 
vers in which the skin buckling might occur before reaching limit load. An 
approach to further investigations of this problem might be to consider an 
approximate strain stiffness matrix of the estimated decrease in skin stiffness 
with buckling, k rr , and formulate 
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> 

-1 





9e ,■ 


8e- 






3Aq r 



- 


3Aq 1 


J 

3Af r 

buckl ed 

J 

3Af r 

' 

[I] + 

Cl] - 

*Vr 


9Af r 
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9Af r 

► 

L_ J 








) 






v. > 


If the k matrix of equation (81) is multiplied by a constant factor, c, 
the equation may be expressed as 





3Aq r 

9Af r 


-1 




(86a) 


That is, the unit matrix is simply replaced by a diagonal matrix of the recip- 
rocal of the change factor c. The acceleration change along k 


K 

% 


= i 





e 6 

r r 



(87] 
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and the strain change along j 


F 1 = 1 - [3Aq k /3A-f r l 
for any flight condition. 

Equation (88) may be displayed as a function of c, the change factor, as 
shown in Figure 14. In a skin change, c would be the proportionality of skin 
thickness of the change. In this manner, the effects of the contemplated 
change at many different fuselage locations in a variety of maneuvers can be 
systematically evaluated to select the most suitable compromise of thickness 
factor. 


rr 


9Aq, 

3Af 


[> 


rr 


E 6 

r r 


L e J 6 Jj 
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Strap, Stringer and Strut Changes . - For a stringer or strap change, as 
shown in Figure 15, the stiffness matrix of equation (71) is diagonal with each 
element being the axial stiffness of the corresponding section of the strap. 
Equation (71) may then be expressed as 



3q k /3Af r 

3e J -/SAf r 


rrJ 


3Aq r 

3Af r 



(89) 


In evaluating the effects of a possible strut addition where there is no 
strut at present, as in Figure 16, there is no existing surface on which to 
place strain gages. A pin-ended tube can be put in the position of the strut 
to activate a differential motion transducer such as a linear differential 
voltage transducer or potentiometer. The term k rr is the axial stiffness of 
the strut in equation (89) which may be written 



_ 

r \ 

\ 


3q k /3Af r 



' — H 



3€j/3A f r j 

^ l/k rr + 3Aq r /3Af r 


( 90 ) 
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STRAIN OR VIBRATION WITH CHANGE 
STRAIN OR VIBRATION WITHOUT CHANGE 



CHANGE FACTOR 


Figure 14. - The change in flight response of any coordinate in 
any maneuver displayed as a function of a change 
factor, such as skin thickness. 



Figure 15. - A strap change of stiffness. 




TECHNIQUES AND PROCEDURES FOR VIBRATION TESTING 
OF THE AH-1G HELICOPTER 


The development of effective techniques for shaking the AH-1G and analy- 
zing the acquired vibration data constituted a major part of the research work 
associated with analytical testing. The validity of the methods employed rests 
heavily upon the consistency between the measured structural mobilities and the 
theoretical models for which these mobilities are derived. This consistency is 
critical since the measured mobilities are used not only for obtaining the 
global modal parameters of the test vehicle, but also to derive mobilities 
which were not measured directly. 

Digital signal analyzers have made it possible to measure the response of 
structures to any physically realizable excitation. However, the interpreta- 
tion of measured structure to specified excitation forces is subject to the 
mathematical model used in the process of analyzing the data. The model may be 
more or less sophisticated, depending on the test data. The dynamic testing of 
a structure like the helicopter poses a number of specific problems. These 
problems are associated with: (1) the size and complexity of structure; 

(2) nonuniform distribution of mass, stiffness, and damping; and (3) the cor- 
rect application of linear vibration theory to the process of data acquisition 
and analysis. 

It has been implied that the techniques adopted for the structural dynamic 
testing are closely related with the theory underlying the vibration analysis. 

A discussion of the specific test procedures must necessarily be preceded by 
a brief summary of the theoretical considerations. This chapter addresses: 

(1) the theory of the generalized linear structure; (2) the principal charac- 
teristics of acceleration mobility data; (3) testing procedures for global 
parameters and the estimation of these parameters; (4) testing procedures for 
obtaining mode shapes and the method of calculating mode shapes; and 
(5) methods for deriving mobilities from modal data. 
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Theory of the Generalized Linear Structure 


The dynamic properties of any structure can always be characterized by a 
relationship between a selected set of motion coordinates and the set of exter- 
nally applied forces, i.e.: 


Motion 


Character of 


Force 

Vector 


Structure 

X 

Vector, 


(91) 


The character of the structure implied in equation (91) will be termed mobility. 
If the motion vector is a vector of displacements/velocities/accelerations, the 
character of the structure is termed displacemerit/velocity/acceleration mobility, 
respectively. 


The central phenomenon of vibration theory is cyclic oscillation, hence 
the quantities that go into equation (91) are generally sought in the frequency 
domain; for example, in acceleration measurements : 

IyU)> = [V (w)]{f (w) } (92) 

where {y(a>)l is the Fourier transform of the accelerations; [Y(w)] is the 
acceleration mobility matrix; and (f ( oj) } is the Fourier transform of the vector 
of generalized forces, compatible with the selected set of coordinates. 

From a measurement standpoint, the jkth element of the matrix E Y ( co) ] 
relates the acceleration measured along the jth coordinate when the only force 
acting on the structure is that applied along the kth coordinate; i.e.: 

Yj (w) = Vj k ( oj) f k (u>) when f . ^ k =0 (93) 


Linear vibration response of a structure may be characterized by the following 
conditions: (1) the response of the structure to random forcing is stationary 

in time (i.e., forced vibrations are steady); (2) the elements of the matrix 
[ Y (oj) ] are functions of frequency only, and depend on neither the motion coordi- 
nates, nor the forcing vector; and (3) the mobility matrix [Y(m)] is symmetric; 
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The foregoing conditions have specific practical implications in vibration 
testing and analysis. The first. condition is necessary for any structure to 
survive continuous operation under arbitrary dynamic excitation. The second 
condition more or less stipulates the type of shake test data that. is adequate 
for analysis based on a linear model of the structure. If the mobility func- 
tions measured for different force levels are not the same, the assumption of 
linearity is not satisfied. This is usually the case when only part of the 
structure may be participating in the response. As the force level is 
increased, more and more of the relevant motion coordinates of the structure 
start to participate in the response. The range of linear response is reached 
only when the measured mobility remains unchanged with changing force levels. 
The third requirement is that of reciprocity. If the shaking and measurement 
stations are interchanged, the same mobility should be recorded, otherwise the 
[Y] matrix will not be symmetric, as required by the linear model. 


It is important to note that, in the foregoing characterization of a linear 
system, no assumptions are made about the nature of the damping mechanisms 
occurring in the structures. All the conditions required for linear modeling 
can be verified in the process of the actual shake test of the structure. 


The relationship between the Fourier transform of the force vector and 
that of the displacement vector of a steadily vibrating undamped multiple 
degree of freedom system can be written as: 


[M] + [K] (y (co) > = fU) 


(94) 


where [M] and [K] are real, symmetric mass and stiffness matrices, respecti vely . 
Thus, the displacement mobility matrix for an undamped system is simply: 


[YU)],. = f- a . 2 [M] + [K] 
u l 


(95) 


The presence of damping in its most general form can be modeled by introducing 
a frequency dependent complex damping matrix into equation (95), i.e.: 


[YU)] d 


O. 2 [M] + [K] + [D R U)] + i^U)] 


(96) 
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It is to be carefully noted that this analytical development has meaning 
only in the frequency domain for the general case of damping. This is mainly 
because the physical quantities that can be used to characterize the arbitrary 
damping of a structure are related to the energy dissipated per cycle of oscil- 
lation. In cases where the time domain, force/motion relationship, repre- 
senting the damping mechanism is known, the damped equations of motion can be 
developed in the time domain and then Fourier transformed into the frequency 
domain. However, taking the inverse Fourier transform of the frequency-domain 
equations that may adequately describe an arbitrarily damped system may not 
yield a time-domain system of equations that makes physical sense. In other 
words, arbitrary damping mechanisms may not be susceptible to a time-domain 
description. Mathematical models, developed from time-domain equations of 
motion, usually fail to identify global characteristics of structures with sig- 
nificant damping. 

D 

In general, the elements of [D (w)] are small, compared to those of the 
[K] matrix. Also, in order for reciprocity conditions to be met and for energy 
to be dissipated, the damping matrix must be symmetric and non-negative 
definite over the entire frequency range. 


For a damped system, then: 


[K] + i[D( M )] 



{y(w) } - {f(w)} 


(97) 


Consider the complex, frequency dependent characteristic value problem: 


[K] +■ i [D(u>)] If) ~ A (w) [M] { (ft } 


(98) 


R T 

where {$} = {<j> } + i{<f> } is the complex characteristic vector which can be 

R .1 

assumed to be frequency independent; a(w) = A (w) + iA (oj) is the frequency 
dependent complex eigenvalue. 


If combinations of |An(w)> 

i J 

tion (98), then: 



and 


/■ \ 

A k (J> {<j,} k exist, which satisfy equa- 
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[K] + i[D(w)] ( 4 >}j 

= Vj(a>){*}J [M] 1 4* ) j 

(99) 

and 



(4«}T [K] + i [ D ( to ) ] { 4> 

= [M]{*} k 

( 100 ) 

7 

{(J)} denotes the transpose of By virtue of the symmetry of the [K], [M], 
and [D ( to ) ] matrices, equations (99) and (100) lead to the following orthogon- 
ality relationships: 

{*}] [M] U> k 

= m j 6 jk 

(101 ) 

and 



fOj f[K] + 1[D(u>)]l {*> k 

= kj + idj (u>)] 6 jk 

( 102 ) 

where : 



m . 

J 

= U)T [H] {.Hj 

(103) 

k j 

= w] [K] 

(104) 

d j U) 

" [D ( cu) ] 

(105) 

6 jk 

= 0 j f k 
1 3 = k 

(106) 


It follows that: 

[K] - [M] + i[D(co)] { 4> } ^ - u^rcij + idj(to) <5^ 

■ (v-> - “ 2 X- « Jk < io7 > 
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If the vectors } j exist, it can easily be verified that only the imaginary 

parts of a.(oj) need be frequency dependent. So that: 

J 


Aj(o>) - Aj + i Aj (to) 


(108) 


Indeed, by post-multiplying the transpose of equation (98) by {$}*, which is 
the complex conjugate of the following equation is obtained: 

J 



[K]+i[D(to)] {<}>}• - A-(to) (<j>}T [M] {<£}. 

. J J J J V. 


(109) 


Similarly, the complex conjugate of equation (98) can be premultiplied by 
{ <f> > A to get: 



[K] -i[0(u)] 


{ *>j - 


A*U) <f)T [M] {*}* 


(110) 


From equations (109) and (110); 


Xj (to) + A*(u) = 2 U}T [K] { 4 ,}j/{«f.)T [M] {<},}* (111) 

and 

AjU) - Xj (to) = 2i {f}T [D(co)]{ ( t,}j/{(t,}T (112) 

The right hand side of equation (111) is a frequency independent quantity. How- 
ever, the right hand side of equation (112) is frequency dependent, establishing 
the validity of the claim made in equation (108). 


A complex L x N modal matrix [$] can be defined such that its jth column 
is the L x 1 vector { <^> > ^ ; j = 1, 2, . ..N, where L is the number of coordinates 
chosen to describe the system and N is the number of modes of the system. In 
principle, N is infinite; in practice, over a given frequency range, only a 
finite number of system modes are necessary. 
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Equation (97) can be rewritten to give. 


■Cy (to) > = [$>] 


0] 


[K] - to [M] + i [D(w)] 


M 


-l 


O] {f (co) } (113) 


and, using the orthogonality relationships, leads to the results. 


{y(w)} = [$] 


1 


f(x" - J) + i A? (to) 
U J J 




[$] (f (co) } 


(114) 


By definition: {y(w)> = [Y(w)]{f (u>) } ; hence. 


N 

CY(to)] = J 

n=l 


(cf) } {(f )} 1 

Y n y n 


1 


(*n - <0 2 ) + i X*(u) 


(115) 


R I 2 

anc * A n^ w ^ have uni "^ s °f (frequency ) and, from physical considerations, both 
R I 

A n and A n (w) are positive. It is, therefore, possible to define: 


,R - 
X n “ % 


(116) 


and 


I 2 

A n (w) = g n (w)^ n 


(117) 


The matrix of modal acceleration coefficients of the nth mode is defined as: 


[A] e — {d, } {A} 1 
J n m n lSM n ^ n 


(118) 


The acceleration mobility matrix and the displacement mobility matrix are 
related by: 


CY(co)] = - a . 2 [YU)] 


(119) 
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Making use of equations (116), (117), (118), and (119), the jkth acceleration 
mobility can be written as: 


N 



M = - 


Z A akn 
n=l 


(1 


w 2 /ft n 2 

+ ig p (w) 


( 120 ) 


In the most general case, the dependence of g (w) on frequency may not be known 
However, it is expedient to take advantage of the fact that the ig n (w) term in 

2 2 

equation (120) is dominant only in the frequency range where ui /fi - 1, i.e., 
near the natural frequency of the nth mode. Thus, any suitable representation 
of g (w) which matches the correct value in the neighborhood of u = may be 
assumed. 


The general form of the jkth element of the acceleration mobility matrix 
can be written as: 

N 2/o 2 

l A.. 2 2 (121) 

n=1 Jk (1 - u3 2 /n 2 n ) + ig n (u) 


V - *jk + "jk * E jk + iE jk 


where + i represents the rigid body acceleration coefficients (E Jk is 

usually very small compared to the rest of the terms in the series and is often 

R I 

neglected); A^ kn = Aj kn + iAj kn is the jkth complex element of the matrix of 

modal accelerations for the nth mode; and g n are the natural frequencies and 
damping coefficients of the nth mode, respecti vely ; to is frequency. 
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Characteristics of Acceleration Mobility Data 


Mo de frequency functions 

The real and imaginary parts of Yj k can be written as: 


Jk 


E R 

jk 


N 

I 

n=l 


A R F R (oi) 
jkn n v ' 


A jkn (u) 


and 


Y?. 

Jk 


E ]k 


N 

l 

n=l 


A jkn F n {u) + A jkn F ^ (u) 


or: 


Y jk 


E jk + iE ]k 


1 A., F (ai) 
L jkn n x ‘ 

n=l 


(122) 


(123) 


(124) 


where the mode frequency functions are defined as: 

2 2 , 2 ,. 2 


:R 


F « M _ ' /n n /n n - 1J 

n M = (, 2 /a 2 - l ) 2 + q 2 
n ' J n 


2 2 

cl , * _ 9 n “ /fi n 

F n - , 2,_ 2 


(“Va n i - l) 2 + g n 2 


and 


R n ^ ~ F n ^ + iF n ^ 


(125) 


(126) 


(127) 


Equations (122), (123), and (124) represent the measured acceleration mobility 

as a linear combination of the mode functions. It is, therefore, important to 

acquire a familiarity with the basic characteristics of the mode functions of 

damped systems and the essential features of their linear combinations. Plots 
"R ”1 

of F (to) and F (to) as functions of frequency ratio for three values of the 
damping coefficient are shown in Figure 17. The polar plots of the complex 
F (to) functions are shown in Figure 18. 
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■^R 


The F (co) function is characterized by two peaks at 


>ln = n n /" + % 2 * ®n Z 1 + 9„ 2 


028 ) 


and 


^ / 1 + g ^ + g / 1 + g ^ ■, 

2n n / - } n a n / y n 


(129) 


while the M (w) function has only one peak at y^l + g 2 


Note that. 


w 2n " w ln 


= 2g n ./ 1 + g. 


(130) 


which increases with increasing damping. 

”R “I 

From the plots in Figure 17, it is seen that linear combinations of F and F 
vary rapidly in the vicinity of the natural frequency, and are either negligible 
or slowly varying with frequency in the regions away from the natural frequency. 


Separated Modes . - Equation (121) carries the basic implication that the 
effects of the structure's modes occurring at different frequencies on the mea- 
sured mobility are additive in the frequency domain. If a mode occurs at a fre- 
quency, in the neighborhood of which the contributions from the other modes of 
the structure are either negligible or are weakly varying with frequency, such 
a mode is said to be well separated. The nature of the measured mobility in 
this frequency range will be dominated by that particular mode. 


Classical Modes . - In the case of a classical mode, i.e., when the system 
mode shape is the same for the damped system as it would be for the undamped 
system the A.^ n is a real number, i.e., - 0, and the real part of the mea- 

sured acceleration mobility will show two turning points for each separated 
mode and the imaginary part will show a single turning point only. For a 
classical mode, equation (130) can be approximated to give an estimate of the 
damping coefficient: g n = (aj 2n - w ln )/n n . 
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Figures 19 and 20 show acceleration mobility measurements obtained from a 
helicopter structure. Two close, but distinguishable modes are present. The 
dominant mode can be seen to be very nearly classical, with double turning 
points in the real, and a single turning point in the imaginary mobilities. 


Complex Modes. - For the general case of none! assi cal or complex modes, 

R I 

both A-. and A., are significant. The measured real and imaginary mobilities 

J Km J Krl j 

of a well separated mode contain linear combinations of both F (oj) and F ( co ) 

in proportions given in equations (122) and (123). In particular, if 
I R 

A^ >> Aj|< n s it is the imaginary part of the acceleration mobility which will 
show two turning points, while the real part will show a single turning point 
only. Figures 21 and 22 show an example of this occurrence in the data measured 
from the AH-1G (the shaking coordinate was vertical at the tail, and the mea- 
surement coordinate was vertical at the nose) between 40 Hz and 50 Hz. 


Coupled Modes . - System modes, occurring in frequency ranges such that 
their mutual contributions to the measured mobility in this frequency range are 
rapidly varying functions of frequency, are said to be coupled. 

Mode Clusters . - A mode cluster (or a cluster of modes) is characterized 
by a group of system modes which are coupled together by virtue of the prox- 
imity of their resonances. Mode clusters are usually separated by regions of 
negligible or slowly varying mobility values in the frequency domain. 

Mode clusters generally have the appearance of single modes in wideband, 
low frequency resolution mobility measurements. Higher resolution data usu- 
ally helps to reveal the modal content of a particular mode cluster. Figure 23 
shows broad band (0 - 200 Hz) mobility of a helicopter vertical tail shake, 
measuring vertical acceleration at the nose. Between mode clusters, measured 
mobility is seen to vary slowly close to the zero value. In fact, what appears 
to be a single mode in the 0 - 10 Hz frequency range is actually a cluster of 
two modes as Figures 19 and 20 (which are higher resolution measurements of the 
same mobility in the 5 - 10 Hz range) show. 
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Figure 19. - Measured acceleration mobility of a helicopter between 5.5 

and 10 Hz. (Shaking vertically at the tail, measuring verti- 
cal acceleration at the nose.) 



' 10.099 

• , . i 

-15.000 m 
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Figure 20. - Data of Figure 19 


plotted 


• 1 i 

MTd 

on the Argand Plane. 
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Figure 21. - Measured acceleration mobility of a helicopter between 38 Hz 
and 52 Hz. (Shaking vertically at the tail, measuring verti 
cal acceleration at the nose.) 



Figure 22. - Data of Figure 21 plotted on the Argand Plane. 



Figure 23. - Measured acceleration mobility of a helicopter between 
2 Hz and 200 Hz. (Shaking vertically at the tail, 
measuring vertical acceleration at the nose.) 

Identification of mode clusters is useful in determining which modes 
should be included when truncating equation (121), since the contributions of 
the remaining modes are either negligible or frequency independent. It also 
helps in identifying frequency segments for higher resolution data acquisition. 

Shake Testing for Global Parameters 

Each elastic mode of the structure is characterized by a natural fre- 
quency and a damping coefficient g n which are global properties of the struc- 
ture. These are the only constants that enter into the mode frequency functions 
They are the same for a given mode, regardless of the response coordinate. The 
first stage of modal testing is to determine the global parameters of the domi- 
nant elastic modes which occur inside the frequency range of interest. 
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The experimental data required for determining the global parameters are 
the continuous frequency plots of a number of mobilities which are considered 
to represent the global vibrational behavior of the structure. For a selected 
set of shaking locations, e.g., tail vertical, tail rotor gearbox lateral, etc., 
the transfer functions between the response coordinate and the shaking coordi- 
nate are measured over the determined frequency range. Typical response coor- 
dinates for such measurements are: (1) nose vertical; (2) wing (right and/or 

left) vertical; (3) center of gravity vertical; (4) tail vertical; and 
(5) horizontal stabilizer vertical . 

The test set up for measuring frequency dependent mobility functions is 
shown in Figure 24. The helicopter is suspended as a free body by soft rubber 
bungee chords. The configuration shown has the shaker located vertically at 
the tail and the response accel erometer at the horizontal stabilizer vertical. 
Signals for driving the electromagnetic shaker originate from the signal gener- 
ator. A force gage, installed at the point of force application generates 
voltage signals which are proportional to the applied force. These signals are 
inputs to the dual channel digital signal analyzer. The accel erometer at the 
response coordinate generates voltage signals, proportional to the response 
acceleration, which are also inputs to the digital signal analyzer. 

The signal analyzer is capable of sampling the time-domain force and 
response signals, digitizing these samples and computing the real-time Fourier 
transforms of the data. It also computes the least squares estimate of the 
frequency-domain transfer function between the input and output spectra, which 
is the mobility between the response and forcing coordinates. All the fre- 
quency functions computed by the analyzer over the specified frequency interval 
can be stored on cassette tapes for future restoration and analysis. The 
oscilloscope allows the monitoring of the time-domain signals emanating from 
the force and response transducers. The frequency counter is used to precisely 
measure the frequencies of harmonic signals when required. 

The accuracy with which the global parameters can be estimated is criti- 
cally dependent on the quality of the data acquired for this purpose. For each 
pair of force and response locations, a random shake is done with the frequency 
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bandwidth set to span twice the frequency range of interest {in this case, 

0 - 200 Hz, since the modes of interest are between 2 - 100 Hz). This is done 
to insure that the modes up to 100 Hz are not coupled to modes occurring beyond 
100 Hz, as may be the case when a local mode is present. For each new shaking 
station, several force levels are tested until the range of applied force is 
reached where the mobility plots do not depend on the force level anymore. 

This is one of the linearity requirements on the mobility plots. Having estab- 
lished the required force level and the absence of local coupling modes at 
higher frequency, another random shake is done, this time with the bandwidth 
set at 2 - 100 Hz. The above procedure is repeated for all the accelerometers 
which have been selected for global parameter testing. 


The ratio of modal acceleration coefficient to damping ( A j| cn /9 n ) varies 

not only from mode to mode, but also from mobility to mobility, for a given 

mode, and the prominence of the various modes of the structure will be different 

in each of the mobilities recorded. That is to say, a given mode i occurring 

at n.j may appear very prominently on mobility Y.^w), while the same mode may 

not be so significant in the mobility Y- (w), where i designates a response 

\J ~ 

coordinate different from k. This will especially be the case if the mode shape 
associated with mode i has a much larger mode element at coordinate k than 
coordinate a. The prominence of a mode may also be due to light damping. Thus, 
by examining the set of broadband mobilities recorded, it is possible to asso- 
ciate each mode i with the mobility where the mode most prominently appears. 


Although it is possible to obtain rough estimates of the natural frequen- 
cies and damping of the structural elastic modes from these broadband mobility 
plots especially when damping is very light (e.g., peaks of the imaginary 
mobility plots, and frequency separation of the peaks in the real mobility 
plots), there are a number of specific considerations why broadband mobility 
data is not suitable for global parameter extraction. Among these considera- 
tions are: 


Measurement Accuracy . - The low frequency resolution associated with broadband 
mobility measurements tends to introduce errors into the measured mobility 
values due to the phenomenon of leakage. Leakage has to do with a spreading 
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of the energy contained at each discrete frequency over a relatively narrow 
band nearby. Although considerable effort is exerted into reducing leakage 
effects (e.g., by appropriately windowing) by the equipment manufacturers, the 
phenomenon still has to be reckoned with when the frequency resolution gets 
below certain limits. For acceptable measurements, bandwidths of about 25% of 
the center frequency have been recommended. 

Parameter Extraction Accuracy . - Also associated with low frequency resolution 
are inaccuracies in the parameter extraction methods due to the frequency 
spacing between successive data points. The polar plot of mobilities, (see 
Figure 18) in the vicinity of a mode, describes a circular arc. Most methods 
for extracting natural frequencies, damping and modal acceleration coefficients 
are based on fitting a continuous circle through measured data and in some 
cases computing the rate of change of the arc length with frequency. Since the 
frequency data is discrete, arcs of the circle are necessarily approximated by 
segments. The error incurred by approximating a circular arc by a straight 
line segment increases as the frequency spacing between successive data points 
increases. Narrow band data, with bandwidth less than 25% of the natural fre- 
quency of a given mode, have been found to yield sufficiently accurate results. 
Initial estimates of the natural frequencies can be obtained from the broadband 
data. 


For sufficient frequency resolution and to minimize leakage, the following 
bandwidths are recommended for use in narrow band testing using the HP5420A 
signal analyzer. 


87 



Natural Frequency, 

^.j (Hz) 

Bandwidth 

(Hz) 

Equal to or greater than 

But less than 

2 

3 

.5 

3 

4 

.781 

4 

6 

1.000 

6 

8 

1.5625 

8 

12 

2.000 

12 

16 

3.125 

16 

25 

4.000 

25 

32 

6.250 

32 

50 

8.000 

50 

64 

12.500 

64 

100 

16.000 


In cases of g greater than .25 use a broader bandwidth. In all cases use the 
natural frequency as the center frequency. 


Swept Sine Testing 

For all the narrow band mobility measurements, the excitation was achieved 
by applying pure sine wave signals to the electromagnetic shaker and varying 
the frequency of the sine waves over the range spanned by the bandwidth. This 
so-called swept sine technique was preferred to other excitation techniques 
over a narrow frequency band. Among other reasons for choosing the swept sine 
technique are that: 

1. The energy input into each measurement frequency is maximum. 

2. By choosing the right sweep speed (see sketch on following page), 
the steady state sinusoidal response of the structure is achieved 
at each measurement frequency. This is one of the assumptions 
made in the derivation of the generalized linear model. 

3. Measurements are more accurate and reproducibl e. 

4. The sampling frequencies and the adequate number of averages 
are more easily determined. 
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5. Good linearity and reciprocity checks were obtained. 

6. High resolution of close modes can be achieved by selecting 
the right sweep speed. 

Consider an undamped single degree of freedom linear system, described by 
the following equation of forced vibrations: 

mx + kx = Fe lwt (131) 

If the forcing frequency coincides with the undamped natural frequency, i.e., 
co , the response of the system is secular, and grows linearly with time. 

Schematical ly : 



or. 


x + rr x 



int 


x (t) = -i 


Ft int 
2mn. e 


(132) 

(133) 


However, because of the various dissipative mechanisms which constitute damping, 
the oscillations reach a limiting amplitude after some characteristic time t. 
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For damped hysteretic damping, g, the steady state response is governed by 


x + (1 + ig)ft 2 x = - e iQt 
m 


( 134 ) 


or. 


x ( t ) = - e int 

mgjr 


( 135 ) 


The steady state amplitude is given by 


Xs = 


mgft 


( 136 ) 


and the characteristic time for reaching steady state response can be estimated 
by equating. 


X T (r) = Xs 

which gives. 



( 137 ) 


(. 138 ) 
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thus: 


0 39) 


_ 2 . 1 
T gn Trgf 

Suppose there are two neighboring structural modes with the natural frequencies 
separated by Af Hz. To resolve these two close modes, the speed at which the 
excitation frequency is changing must be of the order of: 


v = VA f = Trgf A f Hz/sec (140) 

where g = damping coefficient (lower bound) 

f = frequency in Hz 
A.p = mode resolution in Hz. 

If an estimate of the lower bound of the damping coefficient and the required 
mode resolution are available, the sweep speed required for swept sine shake 
testing is directly proportional to the frequency, 

V(A f ) 

i.e. — — = const (141) 


The relationship between the linear scale and the logarithmic scale on the 
signal generator is: 


dec 


= log 


10 f 


hz 

0 


(142) 


where 

f hz = frequency in Hertz; (f Q £ f^ z <_ 10^ ) 

fdec = frequency in decades; (0 £ f dec _< 1.0) 
fg = base frequency on the scale. 


From equation (142) 


f hz ■ f 0 * 10 


dec 


(143) 
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X 


Hz/sec. 


(144) 


df 


bz 


df 


Sweep speed v = ^ 


hz 


df 


dec 


dec 


dec 


dt 


From equation (143) 


df 

df 


hz 

dec 


(fg s.n 10) x 10 dec = f hz dn 10 


(145) 


and 



df 


dec 


dt 


£n 10 


(146) 


Thus, by selecting a constant logarithmic sweep speed (df dec /dt = const. = a), 
equation (141) is automatically satisfied. 


The constant a is determined by substituting the desired value of into 

'hz 

equation (146). For example, if at 2 hz we desire a sweep rate of 0.01 Hz/sec, 
then 


0.01 (60) j ^ , . 

~ 2ln 1 0 dec/mm = .13 dec/min 


(147) 


Estimation of Global Parameters 


Various techniques have been developed for estimating the natural frequen- 
cies and damping coefficients of the elastic modes of a structure from mobility 
data. In all cases, certain assumptions have to be made about these modes. 

The simplest case is when the mode is well separated and lightly damped. For 
such modes, the natural frequency can be approximated by the peak of the imag- 
inary displacement mobility. The damping coefficient can be estimated as: 


2 2 
1 w 2 - W 1 
n = — 


(jl)o — CO n 


(148) 
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where and are the turning point frequencies in the real displacement 
mobility. The above simple case is almost exclusively reserved for simple 
structures with uniform distribution of mass, stiffness, and damping. Very 
few of the modes of the helicopter can be treated this way. 

2 9 

The ds/df method of Kennedy and Pancu . - The following is a more general 
approach which has been found to work well for both classical and complex, 
close or separated modes. By analogy with equations (122) through (127), the 
jkth displacement mobility can be expressed as: 


E. N 
Y,,. = -&■ + l 


jk 2 

-(jj 


n=l 
+ i 


A r f R - A 1 F 1 ' 
[ rt jkn h n rt jkn h n 


«R pi + aI pR 
[ jkn R n jkn R nJ 


(149) 


R I 

where, A^ n = A.^ n + iA.^ n , I s The jkth modal acceleration coefficient of nth 

mode and E-. is the contribution from the rigid body modes. 

J K 

Recall that 


F n (u) 


F F (oj) + i F*(u) £ Fp (“) 


(150) 


which gives 


F n<“> 


1 - io 2 /fi 2 


1 - w 2 /« 2 


+ 9 


(151) 


and 


F 1 ( co ) 
n v 7 


- 9 


n 


1 - 


(o 2 /njj 


+ 9, 


(152) 


In the immediate vicinity of the nth natural frequency, the displacement 
mobility can be approximated by: 
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Y jk ( “ “ R n ) 


A jkn F n<“) 


A jkn + C nV d n 


+ i 


A jkn F ;(“) + A Skn F n^ + C n J2 + d n 

ti 


(153) 


In equation (153), the sum of the contributions from all other modes has been 
represented by a complex straight line: 

““n 


co i iR . • 1 1 

72 + d n + 1 d n 


Dropping the subscripts j, k, n and writing the real and imaginary parts of the 
displacement mobility separately, gives 


y r 

(w - 

- fi) 

= a r 

f r ( w ) 

- A 1 

F^(oj) + 

Y 1 

(u - 

- si) 

. a r 

F J ( M ) 

i — i 
< 

+ 

F R ((jo) + 


R 4 + d R 

(154) 

Si 

1 V d 1 

(155) 


n 


If the nth mode is classical and well separated, the imaginary part of the modal 

acceleration coefficient vanishes and the contributions from other modes are 

I R i 

nearly independent of frequency. In other words, A , C , and C vanish. Thus, 


Y R ( gj — Si) — A R F R (oi) + d R 

(156) 

Y 1 (oj - si) - A R F 1 (to) + d 1 

(157) 


The peak of the imaginary mobility occurs when 

,1 


% - 0 = A R 4 
dw 


or, 


A 


d w 


R 2. /si g 


fRw) 


2 ,^2 


1 - oj' /Si 


jL r 


= 0 


9 9 

1 - w /ST 


.2 


+ g 


(158) 


(159) 
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which is when 


to 2 /Q 2 = 1 , 


as stated previously. 


The peaks of the real displacement mobility occur when, 
,R 


= o = J- 

j 2 U a 2 

dw dai 


F R (uj) 


or. 


(l - uW) 2 - g‘ 


= 0 


1 - w 2 /^ 2 


+ g 


which gives peaks at 


(160) 


(161) 


-i = 1 - 9 


062 ) 


and 


or. 




-4 = i + g 


g = 

as stated previously. 


2 2 

-j 0)2 ~ w-j 0)2 ~ w-| 
- = 

£2 Q 


(163) 


(164) 


When the mode is complex, equations (154) and (155) indicate that both the 

real and imaginary parts of the displacement mobility contain linear combi na- 
R I 2 

tions of F (w), F (oj), oj , and constants. The peaks in the mobilities in the 
general case may not be simply related to the natural frequency and damping 
coefficients. Naturally, different degrees of approximations are feasible, 
depending on how complicated the situation really is. 


A general technique, which has been found applicable to the majority of 
modes encountered on the AH-1G helicopter is based on the rate of change of the 
arc length of the modal curve (plotted on the complex plane, i.e., the plot of 
the Y 1 against Y R with frequency as a parameter). 
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dU 2 ) 


d(<7) 


d(oi 2 ) 


(165) 


where s is the arc length. 

The rate of change of the arc length with respect to the square of fre- 
quency is stationary when 


dY R d 2 Y R + dY 1 d 2 Y r 
d 2 s = n = dU 2 ) d(ai 2 ) d(u 2 ) d(u 2 ) 

df 2 i 2 /? FT 2 F FT 2 

d(“ ) / dY R + dY 1 

/ d(ai 2 )J + d(ui 2 ) 


dY R d 2 Y R dY 1 dV _ n 

2 + 2 u 

d( u 2 ) d(u 2 ) d(w 2 ) d( u 2 ) 


Upon substitution of equations (154) and (155) into equation (167) and simpli- 

d s 

fying, the following condition for the peak of the 2 P^ ot 1S obtained: 

d(co) 


|A 1 2 -V HV 

d(a) 2 ) d(uT) 


'I r r R n H 2 r R T A 1 

dF I , C_ LR _d_F _ A 

2y [ 2 ? 2 A ? 2 

A A d( w 2 ) d(w ) 


pi D A 1 T A R 
, C /i R dr ± /i I dr —n 

+ ~2 A — rr + A — 7i\ - 0 

n d(a> 2 ) d(w) 


? R 2 T 2 

where |A| a (A ) + (A ) . 


For a well separated mode, the constants C R and C 1 will be nearly zero and 
equation (168) gives 
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( 169 ) 


Q. 

f 

r^i 

2 

+ 

r dF 1 1 

CM 

dU 2 ) 


_d(o) 2 )_ 

— i 

CM 

3 

y 


Equation (169) does not involve the modal acceleration coefficients A^ and A 1 . 
Thus, the condition for the peaking of the rate of change of the arc length 
with respect to frequency squared, holds true regardless of how complex the 
mode is, as long as it is well separated. 


Equation (169) can be expanded, making use of equations (151) and (152), 
and the result is: 


d 

dU 2 ) 


2' 

1 - 

1 2 

-g 2 

2 

+ 4g £ 

21 

1 - — 

2 

-l ^ . 



n 


f 2] 

2 

l - 

2 

+ g 2 

l n J 



(170) 


Equation (170) can be evaluated to yield 


d 

d(co 2 ) 


2 


03 


i 

4 L 1 - n 2 J 


f 2) 2 „ 

1 - % + 9 2 

Q 

2 

n 10 

f 2) 

i w j. 2 

i - -o + g 

l n J 

3 


(171) 


ds 

Thus, for a well separated mode, the peak of the plot will occur at the 

dU 2 ) 

natural frequency, regardless of whether the mode is complex or classical. Any 
suitable finite difference scheme can be used to compute ds „ - from measured 


data using equation (165) 


d(w 2 ) 


It turns out that even for modes that are not well separated, the peaks 

of the ~^~2~ Pl°t still give good approximations to the natural frequencies. 
d(o) ) 

To establish why this is so, consider equation (168) term by term. The first 
term vanishes at the natural frequency, as we have already seen. The remaining 
terms can be rearranged as: 
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(172) 


(173) 


Lightly damped modes generally tend to be classical and well separated. This 

is understandable, since in the limit of zero damping, a classical undamped 

8 3 

mode results. Thus, the low damping which will tend to drive 2/fig up, also 
I R R I 

drives C A - C A down, effectively neutralizing the expression. This con- 
sequently reduces the error incurred by approximating the natural frequency by 

n 

the peak of the ds/d(w ) plot. Experience has shown that modes which are too 

o 

close to be resolved by the ds/d(w ) routine may not be resolvable by any other 
method presently known. 


The diameter of the modal circle that fits the curvature of the displace- 
ment mobility plot (on the polar plane) in the vicinity of the natural fre- 
quency is: 



(174) 


At the natural frequency; 

dY R _ - A R C R 

2 2 4 + 2 

d(u ) w = fi g fi fi 


(175) 
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and 


dY' 


d(u 2 ) 




- A 1 

2 4 

g n 


(176) 


Substituting equations (175) and (176) into equation (165) gives 

| A | 

= - a - - v + Error of approximation (177) 

g q 

9 

The plot of arc length s against w has a characteristic S shape, as shown by 
the sketch, below: 


ds 

d( w 2 ) 


S 



By fitting the best straight line to the inflection region of the S curve, one 
obtains a, the gradient of this line. 

a = |A|/g 2 n 4 (178) 

From equations (174) and (178), the damping coefficient is evaluated as: 

g * D/afi 2 (179) 


where D is the diameter of the circle fit to the mobility data, a is the gra- 
dient of the line fit to the S plot and n is the natural frequency. 
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Testing for Orthonormal Modes and Mode Shapes 


The mode shapes of any structure are related to the modal acceleration 
coefficients as shown in equation (118); i.e.. 


A jkn m n ^jn ^kn 


(180) 


where <j> . and <j> 
jn 

nth mode; nr is 
n 


k are mode shape elements at the jth and kth coordinates of the 
the generalized mass of the nth mode. There are two basic 


types of orthonormal modes which can be distinguished by considering the nature 


of the response and the excitation. The ordinary vibration orthonormal mode 
element, v, has units of (1/mass)^ 2 . The products of the jth and kth ortho- 


normal mode elements and the mode frequency function summed over the modes 
defines the jk vibration mobility. On the other hand, the strain orthonormal 
mode element, has units of the square root of the reciprocal of 

force x seconds 2 x length. The products of the jth strain and kth vibration 
orthonormal mode elements and the mode frequency function summed over the 
modes defines the jk strain mobility. The types of orthonormal modes used in 
analytical testing and the corresponding types of mobilities are summarized 


in Table XVIII. 


TABLE XVIII. - SUMMARY OF MOBILITY AND ORTHONORMAL MODE ELEMENTS 
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The orthonormal mode elements are defined as: 


and 


Thus , 


and 


n 


kn r— \n 
✓m 

n 


A jkn ^jn ^kn 


A jjn ^ j n ^ 


It follows from equations (183) and (184) that 

A. 


kn 


jkn 


SK7 


Jjn 


It can also be deduced that 


/ A A 

jkn jikn 

A 0 • 

£jn 


(181) 

(182) 

(183) 

(184) 

(185) 

(186) 


The choice of using either equation (185) or equation (186) to determine 

depends on the accessibility of the modal acceleration coefficients involved. 

Note that two shaking stations are involved in equation (186), whereas only one 

shaking station is involved in equation (185). It may turn out that the 

driving-point data that yields Ajj n are such that accurate estimations of the 

A., for a number of the modes are not easy. This may in part be due to a 
J J ^ 

strong local mode coupling, or a residual effect. In cases where this is so, 
it may be better to shake at a number of coordinates, and then use schemes 
similar to that in equation (186). 

Consistency of the phase angle in equation (.186) is achieved in the fol- 
lowing manner. For an orthonormal mode element, in the nth mode, of large mag- 
nitude, say , let 
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(187) 


kn 


| A jkn 

• 

A £kn 


A £jn 



2 ^jkn + ^£kn ~ ^jn^ 


where $ is phase angle. For any other orthonormal mode element, say p, 


pn 



IT 

A„ 

£pn 

A £jn 



^pkn ” ^kn^ 


(188) 


Mode shapes of the structure, normalized with respect to the highest mode 
element can be obtained directly from the modal acceleration coefficients as: 


= 


V A. 

j ,max,n 


fA ) 

jin 

A. 0 

j2n 


A jNnJ 


(189) 


where A. is the modal acceleration coefficient with the maximum amplitude 

j »max , n 

in the column corresponding to the nth mode, when shaking at the jth coordinate. 
The generalized mass, corresponding to the mode shape thus normalized is com- 
puted from equation (180) as: 


m n ^jn ^kn^ A jkn A jjn^ A j,max,n^ 


(190) 


In order to obtain the elements of the orthonormal modes and mode shapes, 
the modal acceleration coefficients of all the modes for the mobilities 
relating the response coordinates to the shaking coordinates have to be deter- 
mined. The computational scheme for determining the modal acceleration coef- 
ficients requires mobility data at discrete frequencies. The technique, herein 
referred to as the matrix difference method, was developed by F. D. Bartlett, 
Jr., of the Structures Laboratory, USARTL (AVRADCOM). The matrix difference 
method is well suited to processing large numbers of transducers for modal 
analysis using multiplexing data acquisition systems common in the helicopter 
industry. The natural frequencies and modal damping must be determined before- 
hand. 
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I 


For two frequencies cot and co^ in the region of the natural frequency of 
the ith mode, equation (124) could be written thus: 


where 




^jk^ w i ^ ^ A jkn A i F n 


A i F n 






(191) 


Equation (191) can be written for all the remaining modes, having selected the 
corresponding pairs of frequencies. The resulting system of equations is the 
matrix difference equation: 


‘1 



A l' F l 

A 1 F 2 • 

* A 1 F N 


A jkl 

*2 

^k 

>■ = 

a 2 F 1 

A 2 F 2 • 

• A 2 F N 


A jk2 ' 

a n 



a n f i 

A N F 2 * 

• a n f n 


A jkN 

> 


H-k} ■ M 


From which 

{ A j k } - ^r 1 {aV.J 


(192) 


(193) 


(194) 


An immediate observation about the matrix difference scheme is that all 
contributions to the mobilities near a given mode which are weakly varying with 
frequency, such as the effects of distant modes or rigid body modes, are sub- 
tracted out. By proper selection of and a-F can be made such that 
A.j F.j is large and A..F. is small for all jVi . Experience shows that 

<4 ’= fl^l + g i /2) (195) 
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and 

“i = V 1 ' V 2) (196) 

are the most effective choices for the upper and lower discrete frequencies. 

For these discrete frequencies, the matrix [aF] is well conditioned for inver- 
sion, since the off diagonal terms are small compared to the diagonal terms. 

Test procedure . - Figure 25 shows the schematic of the instrumentation 
set-up for the shake test for orthonormal modes and mode shapes. Signals from 
all the accelerometers and from the force gage are transmitted via telemetry 
to a computer where the transfer functions between the response coordinates and 
the force coordinate are computed and printed out. The excitation signals are 
sinusoidal at the discrete frequencies and wT for i = 1, 2, . . . N. The 

same force levels used for the swept sine global parameter shake test are also 
used for the modal shake test at the corresponding discrete frequencies. 

Derivation of Mobilities 

Underlying any technique of modal analysis is the principle of linear 
decomposition of structural response mobility into contributions from the natu- 
ral modes occurring between a chosen frequency interval. The preceding methods 
estimate not only the natural frequencies and damping coefficients of each 
mode, but also the modal acceleration coefficients of each modal contribution 
to the mobility between response and forcing coordinates. 

Subsequent to the determination of the modal parameters and modal con- 
stants, the next logical step is to reconstruct mobilities both between a pair 
of forcing and response coordinates over a continuous frequency interval, and 
at a chosen frequency between several pairs of forcing and response coordinates. 
By comparing the mobility derived over a continuous frequency range with the 
measured mobility over the same frequency range, some assessment of the accu- 
racy of the global parameter estimations can be made. The comparison of dis- 
crete frequency mobilities for a large number of coordinate pairs allows the 
assessment of the acceptabi 1 i ty of the orthonormal mode and mode shape calcula- 
tions. The results of these comparisons build the confidence in the mobilities 
which are derived but not actually measured. 
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from force gage 



Telemetric 

data acquistion system 




Figure 25. - Schematic of test set-up for matrix difference 
method of modal testing. 







Comparison of measured and simulated mo bilities over frequency band. - 
Global parameters and g n of system modes occurring within a specified fre- 
quency range can be satisfactorily estimated using methods based on the prop- 
erties of the mode functions, F (w). The matrix difference method can then be 

R I 

used to calculate the modal acceleration coefficients (A^, A.^) of the rele- 
vant elastic modes. Figure 26 shows plots of the mobility measured between 
0 and 50 Hz. Table XIX summarizes the parameters estimated between 0 and 50 Hz 
from the tail vertical shake/nose vertical acceleration data. Using the 
parameters of Table XIX and equation (121), without including the rigid body 
coefficients, the plots of Figures 27, 28, and 29 were generated. 


The computed and measured mobilities are superimposed in Figures 28 and 29. 
It is seen that the two plots agree to within a frequency independent complex 
constant, which is an estimate of the contribution of the rigid body modes. 


THfMS ft* : 62 M: 1314 [VM 




Figure 26. - Measured acceleration mobility data between 
2 Hz and 50 Hz. 
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TABLE XIX. - ESTIMATED PARAMETERS BETWEEN 0 - 50 Hz 
(TAIL VERTICAL SHAKE, NOSE VERTICAL 
ACCELERATION) 


Mode No. 
n 

Natural 

frequency 

n n 

(Hz) 

Damping 

coefficient 

^n 

Nose/tail modal acceleration 
coefficient 

Real 

a r 

m ZN, ZT, n 
g/1000 N 

Imagi nary 
A 1 

A ZN, ZT, n 
g/1000 N 

1 

7.33 

0.062 

.165 

.032 

2 

8.09 

0.12 

.102 

- .081 

3 

13.3 

0.13 

- .004 

- .001 

4 

15.97 

0.085 

.005 

.032 

5 

16.35 

0.05 

.012 

0 

6 

17.63 

0.08 

.031 

- .013 

7 

22.1 

0.15 

- .017 

- .092 

8 

28.4 

0.11 

- .024 

.026 

9 

40.7 

0.12 

.002 

- .005 

10 

45.3 

0.026 

.009 

.032 


trpns 

Ri: 591 

•fl: 16 





Figure 27. - Numerical simulation of the elastic component 
of the acceleration mobility data. 





Modal Series Method . - Equation (121) can be rewritten in the following 


form: 



+ R >> 


N 



jkn 


F n<»> 


+ Rj k (to) 


(197) 


where 


r!:. (w) is the low frequency mobility residual; i.e., contribu- 
J tions to the mobility by elastic modes which occur at 
frequencies below the lower test frequency limit. 


r H 

jk 


(w) is the high frequency mobility residual. 


The rigid body acceleration coefficient, 
5 

weights data. 


E 


jk * 


is determined from geometry and 


If the lower test frequency limit is near zero, it follows that R L = 0. 

[_[ 

The higher test frequency limit is usually selected high enough so that R can 
be safely assumed to vanish for all but certain driving-point mobilities which 
may suffer either from local mode effects or from high frequency mode residuals. 


When all the global modal parameters (natural frequencies and damping coef- 
ficients) and the modal acceleration coefficients have been determined (see 
Table XX), the acceleration mobilities between pairs of motion coordinates which 
do not include the shaking coordinate can be computed from: 


am 


= E 


m 


n 

+ i 

n=l 


A A 
H £kn mkn 

^kkn 


F n (u) 


(198) 


where k is the coordinate of the shaking station for the data which generated 
A ^ and A It is necessary to select the shaking station k such that there 
is no local mode or high frequency residual effect on the estimated value of 

\kn* 


If only of the modes are well defined by shaking at k, while the 
remaining Np modes are better defined by the shake at p, then 

IV a A A A 

;; _ p- , r £kn mkn p , \ . r K xpn mpn p , x 

+ Jl F " W + v F " (ui) 


(199) 
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TABLE XX. - SUMMARY OF ESTIMATED MODAL PARAMETERS 
FOR AH-1G HELICOPTER. 



no 






















































SPECIAL CONSIDERATIONS IN MODAL ANALYSIS 


Shaking Locations 

Nonlinearities such as gradually hardening stiffness and ankylotic stiction 
are alleviated by using sufficient shaking force so that the mobility is inde- 
pendent of the shaking force. The power input from the shaker is equal to the 
power dissipated by the helicopter through damping and is proportional to the 
product of the square of the shaking force and the imaginary velocity mobility 
(or real acceleration mobility divided by frequency). Maximum motion for a 
given shaking force occurs at shaking point locations of high mobility. These 
are generally at structural terminations such as high on a fin, at wing tips 
or the hub. 


It is also important to shake in the vicinity of any possible local modes 

of interest. Consider the case in which, for a given mode, ^ = .2 and ^ = 2.0. 

For shaking at j, A.. = .04 and A.. = .40. The response near this natural fre- 
J J KJ 

quency may be unnoticeable at the driving point, j, but evident in the kj trans- 
fer mobility which indicates that k, or a point near it, should be a shaking 
point. In this example, A ^ is two orders of magnitude larger than A... 


Once modal accelerations have been obtained for and among the stations at 
which the aircraft is excited, one shaker at a time, all the modal accelerations 
can be used simultaneously for the best estimate of the orthonormal mode ele- 
ments. Although the number of shaking locations depends on many factors, it 
appears that five or six should be sufficient for most helicopter uses. For 
the jth motion coordinate and N shaker locations the absolute value 


A jl + l A j2 1 + 


+ A jN 


/A 


11 


72 


+ . . . + 2 A-j ^ + 2A 13 


.. + 2 A 23 + 


A 



N 


I 

i=l 
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The phase angle of the jth orthonormal mode element is given by 
1 N 1 N 

= N J, +ji - 2N 1 | 1 +ii (201) 

2 

If the ds/df method is used to obtain the modal accelerations at and 
among the shaking locations and the matrix difference method for all modal 
accelerations, then 

N 

I |A 11 /g| 
i=l I J' I 

— | (202) 

A,,/ 9 2 + A 22 /g 2 + ... + 2 |A-, 2 /g 2 + 2|A 13 /g 2 | + ,.T 

High Frequency Residuals 

The authors discovered that high frequency (e.g., 400 Hz) local modes may 
cause distorted translational mobility plots, with force shaking at low fre- 
quencies (e.g., 10 Hz) at or near a driving-point on the AH-1G. Dr. Albert 
Klosterman and Dr. Jason Lemon of Structural Dynamics Research Corporation 
noted high frequency residual effects and found the means of correction in 
testing a 1750 HP electric motor for the United States Steel Corporation. 

Dr. David Ewins of the Imperial College of Science and Technology, London, 
found the same cause and effect in rotational mobilities measured on a turbine 
rotor. ^ The effect does not occur at all driving-points, but it is not uncom- 
mon. It occurred in tail shaking and vertical shaking at FS Z400 on the AH-1G 
as shown in Figure 30. It did not occur in hub shaking or shaking at the tail 
rotor gear box on the AH-1G. 

The effect is associated with shaking points of very small mass concentra- 
tion and elastic extension as in the cantilever beam of Figure 31. In a heli- 
copter these correspond to soft, but not weak points, as opposed to hard points 
usually, but not necessarily, such as skin attachment points or tail skids. 

In some cases of helicopter fuselage concern, the high frequency local 
modes are inconsequential in engineering and should be ignored in modal 
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Figure 30. - Mobility curves with the high frequency residual effect 
measured on the AH-1G. 



superposition for mobilities. If the point in question is one at which there 
is a contemplated external stores change, the point would usually have to be 
hardened for the attachment and the effect would disappear. Only when true 
very local response is desired should the high frequency local modes be 
retained and this would be the case on panels and decks. If the high frequency 
residual effect occurs in direct shaking at a change point which is to be har- 
dened and the effect not removed then there could be significant errors in the 
analysis of the effects of the change. 
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Figure 32 shows a simple chain having the residual effect. The natural 
frequencies are at 4.95 Hz and 123.28 Hz, the latter being a local mode of the 
small mass. The acceleration mobility from 1 to 10 Hz for the response at 1 
due to force at 1, Y(l,l) shown in Figure 33 shows the residual effect. 



969 N/m 300 N/m 


Figure 32. - Simple chain system with 5 % structural damping. 



Figure 33. - Driving-point acceleration mobility at mass 1 
of the chain of Figure 32. 

Y(1 ,1 ) is almost identical to Y(l,2) and would be indistinguishable from 
Y(2, 2) , shown in Figure 34, if the local mode were removed. The slope of the 


114 





mobility from 4.8 to 5.07 Hz is virtually the same in Figure 33 as in Figure 34 
and differential mobility parameter extraction methods, such as dS/df^ and the 
matrix difference method, on the 4.95 Hz mode would be negligibly affected by 
the residual effect. These methods can be used to accurately obtain modal 
accelerations at driving points whenever the slope (dY/dw) resulting from high 
frequency local modes is small compared to the maximum real or imaginary mobil- 
ity slope at the natural frequencies of interest. The accuracy of parameter 
extraction methods that are not of the differential mobility type is seriously 
affected at driving points exhibiting the high frequency residual effect. 



Frequency, Hz 


Figure 34. - Driving-point mobility of mass 2 of the 
chain of Figure 32. 


If the high frequency residual effect is to be retained there are two ways 
to deal with it without shaking to very high frequencies. Because it is essen- 
tially a driving-point phenomenon, one can shake at each of the driving points 
having the effect and represent the effect as a flexibility term equal to the 
constant difference between measured displacement mobility and the displacement 
mobility synthesized over the frequency range of testing as shown by Klosterman 
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and Lemon^, and Ewins 11 . Alternatively, one could shake at only one point with 
the points having the effect, and at which driving point mobilities are desired, 
weighted by attached masses. The effect of the attached masses could then be 
removed using the analytical testing equations as follows: 


11 


22 


33 


[I] - 


[Y rr ] 




[m rr ] 


-1 


Y (s) 

Y n 

Y (s) 

'22 

Y (s) 

33 


( 203 ) 


where is the matrix of measured acceleration mobilities at and among the 
motion coordinates with the attached masses, m is the diagonal matrix of the 
attached masses (in three spatial directions if necessary) and the parentheti- 
cal superscript s meaning derived from modal synthesis over the frequency range 
of testing. 


Effect of Damping Estimate Variations in the 
Matrix Difference Method 

The matrix difference method involves the difference between mobilities 
above and below each natural frequency. Let u | be the frequency below the 
xth mode and w £ be the frequency above the xth mode. For illustration pur- 
poses let 


a>xi = and 


w x2 ~ ^ + s 


where s is much less than unity. Then 


Am 


Atu x 5 w xl " w x2 


= (/Ts - fi+s) a = - so or - s *— 


and in terms of displacement mobility: 


( 204 ) 


( 205 ) 
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Y jk^ w xl^ " Y jk^xz^ 


(206) 
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(208) 


In view of equation (205), the xm th term in the gAy matrix may be expressed 
as 


- 2g 


Aw x S 


g M Ay xM 


M ft 2 
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- 9m + 12g, 


1 - 


(209) 


In the x row of the gAy matrix this term shows the coupling of the xth 
mode with the mth mode when ratioed to the diagonal term g Ay. It will be 

X XX 

shown that: (1) when the modes are widely separated (i.e., for any mode x, 

- 1 << (1 - g x ) ft x and >> (1 + g x _-| ) ft x _-, and ft x << (1 - g x+] ) ft x+1 and 
ft + 1 >> (1 + g ) ft ), then a frequency difference of Aw x /ft x - 9 X is preferred 
in that it is least sensitive to errors in g; (2) for close modes, Au> x /ft x should 
be less than g : and (3) for Aw /ft much less than g v , the maximum error in the 
diagonal term will be proportional to the error in g . 

Figure 35 shows a plot of the x row real part of the gAy matrix for hys- 
teretic damping of 10% for various values of Aw x /ft x equal to or less than 
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g x - .1. As seen from Figure 36, the imaginary part of the coupling of the x 
mode with all other modes also drops off with decreasing Aw /fi . This also 
occurs in Figure 35. Figures 37 and 38 illustrate the severe coupling effects 
when Aoo x /^ x is greater than g by a large amount. Note the development of two 
peaks when Aw x /^ x is .4 for g = .1 in Figure 37. 

The matrix difference method remains exact for values of Aw x /ft x much 
larger than g but the increased coupling with nearby modes tends to make the 
gAy matrix less well conditioned both physically and numerically. 

As an illustration of the variation in coupling of two modes in the g 

2 2 

matrix with variation in Aw x /^ x , consider a case of g = .1 and <^ x = (1 - g^)^ . 

The following calculations are made for = 1.05409 and g = .1. 


4 “x /fl x 
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Figure 37. - The real part of the row for the natural frequency in the gAy matrix 
when Aw x /ft x > g x . 
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If g were .01 the same case of close modes has the following coupling 
terms. 
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Sensitivity to g 


When the diagonal term from equation (209) reduces to 


x 


g Ay 
J X M XX 


2g x A» x /n x 


( 210 ) 


AwX 


+ 


As seen in Figure 39, the diagonal term in gAy is least sensitive to the value 
of g when Aw x /ft x - g. In that vicinity a 4:1 variation in g produces only a 
+ 10% change in the value of the diagonal term. 


If Aw x /^ x is very small in equation (210) compared to g then 


Ay 


XX 


+ 2 


AM x /n x 

9x 


( 211 ) 


and percentage errors in the diagonal term become proportional to percentage 
errors in g. 
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If g is so small that there is negligible damping coupling with other 
modes then one can assume any small value for g. 

A(0 x 

Let - — = .05 as an illustration of a case with small g. Then, 
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Ay 

M xx 

Ay xM for S! m 
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0 ° 
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0 . 7 ° 

.005 
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38.46 
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.02 

34.48 
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Effects of Close Modes 

Modes i and i+1 are considered to be close when 

fi i+ l < (1 + 9^ ^ or n. > (1 - g i+1 ) n i+1 (212) 

where g is structural damping and the natural frequency n. + .| > n. . 

It should first be noted that ft. and g. are global properties. Therefore, 
one might find ft. and g^ from one mobility plot and ft^ + ^ and g^ + ^ from a dif- 
ferent mobility. One must be very cautious in dealing with what are presumably 
close modes and make certain that they are, indeed, different modes if the 
modal information is used to synthesize mobilities at which the structure was 
not shaken. The method for doing that is beyond the scope of this report. 
Structures are not purely linear systems and slight differences in measured 
natural frequencies from mobility to mobility do not necessarily indicate dif- 
ferent modes. 
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The jkth displacement mobility of a free system at frequency w is given 


by 


E.. N A... 

V“> ----T - + X -fr 

W 1-1 Si- 
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where Is the jkth rigid body acceleration coefficient. At 1 Hertz, 

Ejk/^4 = Ej k /1559, so for practical purposes in such structures as helicopters 
and airplanes 
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As shown by Kennedy and Pancu' 
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Let 


^jk 1 


dr 


J 1 


dS ik 

be the value of — sA at the ith natural frequency. Then 
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and equation (218) can be written as 
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a limiting case of the matrix difference method. It is not necessary for any i 

d s 1k 

that — k— be a peak. If the modes are not very close, as defined by equa- 
df 2 2 

tion (212), then the g y matrix approaches the unit matrix. In any case, the 

2 2 

diagonal terms of the real g y matrix will always be -1.00 and all other terms 
will have an absolute value less than unity. 


Defining 
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The orthonormal modal elements are obtained in the usual manner and any 
mobility 


Y Mn (u) = - -f 1+ J. ¥ n' 1 m 9 i w i 
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or, in terms of acceleration mobility 
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2 2 

The advantage in using g y matrix for a system with the following natural 
frequencies and modal structural damping coefficients is shown below. 


Damping g 

.1 

.2 

.1 

.15 

.1 


Natural frequency (Hz) 

5 

9 

10 

11 

15 


Forcing frequency 
(Hz) 

5 

- 1 .000 

.065 

.017 

.032 

.012 


9 

.002 

- 1 .000 

.123 

.112 



10 

.001 

.067 

- 1 .000 

.062 

.029 

REAL 

11 

0 

.100 

.117 

- 1 .000 

.041 

2 2 
g u 

15 

0 

.012 

.006 

.028 

-1 .000 


5 

0 

- .041 

- .005 

- .013 

- .003 


9 

0 

0 

- .179 

- .128 

- .007 


10 

0 

.416 

0 

- .423 

- .011 

IMAGINARY 
2 2 

11 

0 

.098 

.144 

0 

- .018 

g y 

15 

0 

.003 

.001 

.010 

0 


5 

1 .000 

.077 

.017 

.035 

.013 


9 

.002 

1.000 

.217 

.170 

.024 

MAGNITUDE 

10 

.001 

.421 

1 .000 

.428 

.031 

2 2 
g u 

11 

0 

.141 

.186 

1.000 

.045 

15 

0 

.013 

.006 

.030 

1 .000 



The 5 Hz mode is seen to be essentially decoupled from the other modes in 
gV. That is, it is decoupled insofar as dS/df 2 is concerned. The 9 Hz mode 
is influential at 10 Hz since 9 Hz x (1 + .2) = 10.8. In the magnitude table 
of g 2 y 2 the 9 Hz mode has a value of .421 at 10 Hz. The coupling, it is seen, 
occurs in the imaginary g 2 y 2 . Note that the 10 Hz mode is only half as influ- 
ential at 9 Hz as the 9 Hz mode is at 10 Hz. Similarly, the 11 Hz mode is 
2.3 times as influential at 10 Hz than the 10 Hz mode at 1 1 Hz as equation (212) 
indicates. In other words, 11 x (1 - .15) = 9.35, less than 10, but 
10 x (1 + .1 ) = 11 , just equal to the next higher mode. 
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2 2 

The g u matrix is clearly well conditioned. If the 9 Hz mode were, 
instead, 9.9 Hz the .421 value in the magnitude table would become .999 at 10 Hz 
and the value at 9.9 Hz from the 10 Hz mode would be .962 and the coupling 
would appear mainly in the real matrix. The conditioning is worsened but not 
necessarily unacceptable. 
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CONCLUDING REMARKS 


Analytical testing provides a practical methodology for combining struc- 
tural dynamic analyses with flight vibration performance. This report has pre- 
sented analytical methods for evaluating the flight vibration and strain effects 
of structural changes. The changed flight responses are characterized by the 
dynamics of the change, original flight test data, and shake test data. Basic 
analytical testing equations were derived to accommodate a category of changes 
which included mass, stiffness, absorbers, and active vibration suppressors. 

Illustrations of the analytical testing methodology were presented using 
flight test and shake test data measured on an AH-1G helicopter. The applica- 
tions were directed to the practical acquisition of helicopter vibration data 
and the utilization of the method. Results of this investigation indicate the 
possibility of combining flight data and shake test data with structural 
changes in the helicopter airframe to predict the effects on flight vibrations 
and strains. The results do not represent experimental verification of analy- 
tical testing in any particular application. 

Modal analysis based upon the measurement and past test processing of 
transfer functions constituted a major part of the research associated with 
analytical testing. The vibration testing procedures and modal analysis tech- 
niques of the AH-1G were described. In addition, special considerations were 
discussed to establish consistent data acquistion and to minimize problems 
associated with modal analysis of complex structures. 


130 



APPENDIX A 


DERIVATION OF THE BASIC ANALYTICAL TESTING EQUATION 


Consider two structures with degrees of freedom F and C connected at I 
interface coordinates. 



Forces and 
momen ts 


The structure F is excited by several external forces and moments, referred to 
as forces for brevity. Let {q} define the response of structure F alone such 
that 


Z FF 

Z Fl’ 

« 

q F 

. = 

> 

f F 

Z IF 

M 

i — i 
i — i 


*1 


0 

_ 

j 


, 


, > 


(Al) 


The structure C represents a dynamic change and the coupled equations of motion 
are, assuming the external forces do not change 


Z FF 

Z FI 

0 


q F 


> 

f F 

Z IF 

( Z II + Z lP 

L IC 

■ 

q i 

. — ■ 

0 

0 

L Cl 

z^ 

^cc 


q c 


0 

. J 


where the superscript C defines the free body impedances for structure C and 
{ q ' } defines the response of the coupled system. From the definition of the 
matrix inverse it follows that 
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0 


1+1^ 

L n L u L u 


u ^ci ^ccj |JCF 

Equation A3 gives the following relationships 

Z rr Y ' r- + Z. T Y' = T I 1 


FF 

y ff 

+ 

t-H 

U_ 

M 

-< 
1 — 1 - 
T| 

IF 

y ff 

+ 

< Z II 

+ Z II 

C 

Cl 

Y iF 

+ 

z c 

cc 

-< 
o - 
~n 


Substituting for Y£ F from equation (A6) into equation (A5) and simplifying 
leads to the result 


Z IF Y FF + < Z II + Y II ) Y IF = 0 


c 7 C 7 C-1 7 c 
L U " IC CC ^CI 


Equation (A8) is obtained by considering the matrix inverse definition for 
structure C alone such that 


k 1 
I 

\ 


Solving for Yjj from equation (A9) gives equation (A8). Considering equa- 
tion (A3), it also follows that 
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Z FF 

i — i 

- u_ 
>- 

o 

II 

1— 1 
- 1 — 1 
>- 

1 — 1 
u. 

M 

+ 

(A10) 

Z IF 

y fi 

+ ( z ii + z ii) Y h + z ic Y ci f 1 j 

(All) 

L Cl 

Y h 

+ Z CC Y CI ■ 0 

(A12) 


Substituting for Y£j from equation (A12) into equation (All) and simplifying 
leads to the result 

(A13) 


Y 1 

IF Y FI 


(Z 


II 


+ y?; 1 ) y; 


ii 


ii 


I 


Combining equations (A4), (A7), (A10), and (A13) gives 


i — i 

Ll_ 

M 

LI- 

LL- 

IE- 


y ff 

Y 1 
Y FI 


s: 

Z Z + Y^ - ^ 

L ll r II 


Y iF 

l — i 

— » — i 
>- 

ii 

I 


or, after rearranging equation (A14) 


" y ff 

Y' 

y fi 


y ff 

y fi 


1 — 1 

1 

o 

% 

Y h 


y if 

Y II 


Y C-1 Y 
Y II y if 

L. 

C-l 

I + Y Y 

1 Y II Y II 


(A14) 


(A! 5 ) 


Making use of the matrix inverse definition, matrix algebra reduces equa- 
tion (A15) to 


" y ff 

y fi’ 


y ff 

y fi’ 


1 — \ 

1 

1 

o 

Y 1 

y if 

— 

Y h 


y if 

Y II 


"( y ii + Yii)' 1 y if 

t— 1 
+ 

-< 

i— i o 

HH 1 | 

1 

-< 

1 — 1 
1 — 1 

1 

1 


(A16) 
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1 

FF 

L/U 

►— i 
>- 

1 

i — i 

I T 

>- 

+ 

1 1 

O t-H 
>- 

1 — 1 
Ll_ 
>- 

LX. 

U- 

>- 

II 

(A17) 

IF 

y if “ Y II ^ Y II + Y II^ 1 y if 

(A18) 


Since the external forces are on structure F, postmul tiplying equations (A17) 
and (A18) by the force vector yields 



or, in a slightly different form 



(A20) 


Equations (A 7 9 ) and (A20) are independent of the number of degrees of freedom 
of structures F and C. In other words, the change in response at any point on 
structure F only requires information at that point and at the interface points. 
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APPENDIX B 


COUPLED ROTOR/ FUSELAGE VIBRATIONS AND LOADS 
Consider the helicopter dynamics representation shown below. 



Let z' , z" , and z define the equations of motion (impedances) of the rotor 
alone, fuselage alone, and coupled rotor/fuselage, respectively. Similarly, 
let y 1 , y" , and y define the corresponding rotor alone, fuselage alone, and 
coupled rotor/fuselage responses. The matrix equations of motion can be 
expressed in terms of impedances as: 


Rotor alone 
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I 


Coupled rotor/fuselage 


Z RR Z RH 0 


V 


ff 1 

t r 

Z HR Z HH Z HF 

- 

y H 

= 

0 

0 Z FH Z FF 


y F 


0 


Z RR z rh 0 


y R ' 


f f 1 

r R 

Z HR ^ Z HH +Z HH^ Z HF 

- 

y H 

It 

” 

0 7" 7" 

u z pH Z FF 


y F 

J 


o 

V 


Equation 


(B3a) can be rearranged in the form 


(B3a) 


(B3b) 


Define 


Z RR 

y R + 

Z RH 

y H 


f R 

Z HH 

Z HF 

. 

y H 

... 

Z HR y R 

> 

Z FH 

Z FF 


y F 


0 0 


f H " Z HR y R 


(B4a) 


(B4b) 


(B5) 


where f^ are the rotor hub loads applied to the coupled rotor/fuselage to give 
the coupled hub and fuselage responses as shown below. 
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By definition from equation (B4b) 


f H = Y HH 


From the definition of mobility 


Z HH Z HF 


ry y i 

t hh t hf 

II 

3 

o 

f 

Z FH Z FF 


Y Y 

t FH t FF 


0 I 


Solving equation (B7) shows that 

Y HH = ^ Z HH " Z HF Z FF Z FhP 
Thus, equation (B6) becomes 

f H = ^ Z HH ‘ Z HF Z FF Z FH^ y H 


Since z^ = + z^, equation (B9) can be written as 


f — ( 7 1 4 - 7 11 — 7 11 7 " ^ 7 11 1 V 

t H lz HH Z HH Z HF FF FH ; y H 


From the definition of mobility for the fuselage alone 


z 11 z" 

Z HH Z HF 


”y 1 1 y n 

t hh t hf 


CD 

Z M 7 n 

Z FH Z FF 


yii yn 

t fh t ff 


0 I 


Solving equation (BIT) shows that 


Y" = ( 7 m - 7 " 7 11 - 7 n ) 

Y HH lZ HH Z HF Z FF Z FH J 


and equation (BIO) reduces to 

f H = ^ Z HH + Y HH h y H 


(B6) 


(B7 ) 


(B8) 


(B9) 


(BIO) 


(Bll) 


(B12) 


(B1 3) 
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Substituting for y R from equation (B4a) into equation (B5) also gives 

f H = " Z HR Z RR f R + Z HR Z RR Z RH y H ^ B14 ^ 

Equation (B14) can be simplified by considering a rotor to a fixed hub so that 
equation (Bl) becomes 


Rotor to fixed hub 


7 1 

Z RR 

Z RH 

- 

( * 
y R 

= ■ 

f R 

7 1 

Z HR 

Z HH 


0 


f H 


Solving for f^ from equation (B15) gives 

■f 1 = 7 1 w 1 = 7 1 7 1 ^ -f- 

t H Z HR y R Z HR Z RR T R 

The forces transmitted to a fixed hub are 

f = - f 1 = _ 7 1 7 1 ^ f 

T 0 t H z HR z RR R 

Since z RR = z HR and z RR = z RR , equation ( Bl 4 ) becomes 
f H = f 0 + Z HR Z RR Z RH y H 


(B15) 


(B16) 


(Bl 7) 


(B18) 


Some additional observations concerning fg can be made by considering equa- 
tion (B3b) rearranged as follows: 


Z RR Z RH 


y R ' 


0 


ff 1 
t R 


■ 


> + 


n 

u_ 

>> 


Z HR ^ Z HH +Z HH ) 


y H 


z" 

Z HF 


0 

_ 


< J 


- 


.. > 


(B19a) 


0 




y H 


( z p F ) 


0 


(Bl 9b) 
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Substituting for y F from equation (B19b) into equation (B19a) gives 


■RR 


•RH 


7 1 (z‘ + 7 11 — 7 n 7 11 7 m i 

Z HR vz HH Z HH Z HF Z FF Z FH J 


Further simplification of equation (B20) results in 



- 

y R 


f R 

- 


II 





0 

J 



. J 


^ Z HH ' Z HR Z RR 1 Z RH^ + ^ Z HH " Z HH Z FF^ Z FH^ y H Z HR Z RR ] f R 


(B20) 


(B21) 


Equation (B21 ) can be reduced still further by taking advantage of the mobility 
definition for the rotor alone. 


(B22) 


7 1 7 1 

Z RR Z RH 


Y i v « 

t RR t RH 

hi 

I 0 

Z HR Z HH 


Y ' Y 1 

y hr t hh 


! 1 

°_l 


Solving equation (B22) shows that 


Y 1 = ( 7 1 _ 7 1 7 1 ^ 7 1 ) 

t HH {Z HH Z HR Z RR z RH' 


f 0 ^HH 1 + Y HH ] ) y H 


(B23) 


Substituting equations (B12) and (B23) into equation (B21 ) gives 

^HH 1 + Y HH ^ y H = Z HR Z RR f R ^ B24 ^ 

and equation (B17) can also be written as 


(B25) 


Substituting equation (B25) into equation (B18) and making use of equation (B23) 
results in 


f H ^ Z HH + Y HH ^ y H 


(B26) 
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which is identical to the result derived in equation (B13). Eliminating 
from equation (Bl) gives 


V .-l w* =_7 l 7* ^ f 

Y HH y H Z HR Z RR R 

so that 

f 0 = Y HH ] y H 

Combining equations (B25) and (B28) shows that 

( Y HH^ + Y HH ^ y H = Y HH y H 
or, 

Y HH y H = Y HH ^ y H " y H^ 

Therefore, f H can be expressed as 

f H = f 0 + Z HR Z RR Z RH y H 
f H = f 0 + ^ Z HH " Y HH ^ y H 
f H = (z HH + Y HH^ y H 
f H = Z HH y H + Y HH ^ y H " y H^ 
f H = Y HH ] y H 


(B27 ) 
(B28) 

(B29a) 

(B29b) 

(B30a) 
(B30b) 
( B 30c ) 
(B30d) 
(B30e) 


where 


,,-In-I 


Y HH " ^ Z HH + Y HH ^ 


(B31 ) 


In summary, the rotor hub loads, f H> can be expressed in terms of loads trans 
mitted to a fixed hub and modified by the rotor hub impedances alone, as well 
as the coupled rotor/fuselage hub responses. 
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From another point of view equation (B3b) can be written as 


Z RR + Z RH *H = f R 


7 11 

Z HH 

7 11 

Z HF 

• 


- “ » 

- Z^ R y R - Zj^ H y H 

7 11 

|_ Z FH 

Ll_ 
= LL_ 
N 


y F. 


0 


(B32a) 


(B32b) 


Define f R = - z^ R y R - z‘ H y H (B33) 

where f^ are the rotor hub loads applied to the fuselage alone to give the cou- 
pled hub and fuselage responses. 



Comparing equations (B5) and (B33) shows that 
= f H " Z HH y H 

Applying equations (B30), f^ can be expressed as 
= f 0 ' Y HH 1 y H 

f = Y" ”1 y 

t H y HH y H 

= Y HH ] ^ y H " y H^ 

Also from equation (B32b) 


n= 

>> 

II 

Fyll ' 

t hh 



Y" 

_ p H 


(B34) 


(B35a ) 
(B35b) 
(B35c) 


(B36) 
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Equations (B35) and (B36) serve to explain the nature of coupled rotor/fuselage 
vibrations with respect to rotor/fuselage dynamic interactions. In particular, 
the impact of fuselage dynamics on f^ and consequently on fuselage vibrations. 
Helicopter rotor excitations are very large and one might suspect that hub 
responses of a rotor with no fuselage (y^) may be very large when compared with 
hub responses of the coupled rotor/fuselage (y H ). In other words, 

Yh » Y h (B37) 

As a consequence of equations (B37) and (B28), equation (B35c) becomes 

f H s Y HH y H " f 0 ^ B38 ^ 

Therefore, the rotor hub loads acting on a fuselage without the rotor are nearly 
independent of the fuselage dynamics. Under this restriction the rotor hub 
loads are equal to the forces transmitted to a rigid support. Consider the 
extreme case of a fuselage with zero damping which has a natural frequency pre- 
cisely at a harmonic of the blade passage frequency. Then y^ = y^ since the 
fuselage has zero hub impedance (or infinite hub mobility) in this extreme case. 
If the undamped fuselage has a natural frequency only slightly removed from the 
excitation frequency, then y^ >> y^ appears plausible and the conclusion of 
equation (B38) follows. Also, if with fuselage dynamic changes y^ remains 
small compared to y^, then there is little change in the transmitted forces as 
concluded in equation (B38). 


However, the condition that y^ is small compared to y^ is a sufficient but not 
necessary condition for small changes in f^. From equation (B20) 


r *\ 


Z RR Z RH 

-1 

V 




■ 


>. J 


Z HR ^ Z HH + Y HH^ 


0 

N. J 


(B39) 


Consider the case of a fuselage with a natural frequency at the blade passage 
frequency for a system without the rotor. In this example YJJm is very large 
or Yj^ is very small and equation (B39) indicates that 
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y H = *H 


(840) 


and from equation (B35b) 

f H =Y“- 1 y H (B41) 

If Yj^ does not change appreciably due to fuselage dynamic changes and remains 
small, then there is negligible change in y^ as equation (B39) shows. As a 
result there is negligible change in the transmitted rotor loads, f^, given by 
equation (B41). To illustrate this case, consider an isolated hub-transmission 
where 


Y 


HH 


>> 


Y" 

Y HF 


(B42 ) 


and Yj^ is very small. For a change in the Fth degree of freedom of the fuse- 
lage due to an impedance change AZpp, the changed driving-point hub mobility 
is given by 


aY" = V" 
at HH y HH 


vii V” 

y hf y fh 


AZ 


FF 


+ Y" 
Y FF 


(B43) 


Since Yj^ >> YJ^p, it follows that 


AY" - Y" 
ay HH 


(B44 ) 


and aY^ 1 


remains 


very smal 1 . 


Now consider the opposite situation where Yj^ 
small compared to the z' matrix in equation (B39). 
compared to y^ and the conclusion in equation (B38) 


is very large or Yj^ is very 
This makes y H very small 
f ol 1 ows . 


If at blade passage harmonics, equation (B42) does not hold for various 
fuselage degrees of freedom, then the fuselage has a major as opposed to local 
resonance and the transmitted forces may change considerably due to fuselage 
dynamic changes. However, if the fuselage were separated from the rotor, then 
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such a condition could not be tolerated in practice. A major resonance has 
modal accelerations relative to points other than the hub which are larger 
relative to modal accelerations of other modes. A large modal acceleration 
can be looked upon as a small effective mass since 


^jn ^kn 
jkn “ m n 


(B45) 


Physically, a 
cantly than a 


small effective mass 
large effective mass 


contributes to the response more signifi- 
for the same loading conditions. 


144 



REFERENCES 


1. Schrage, D. P.; and Peskar, R. E.: Helicopter Vibration Requirements. 

Preprint No. 77.33-33, American Helicopter Society, 33rd Annual National 
Forum, May 1977. 

2. Veca, Angelo C.: Vibration Effects on Helicopter Reliability and 

Maintainability. USAAMRDL TR73-11 , April 1973. 

3. Berman, A.: A Generalized Coupling Technique for the Dynamic Analysis of 

Structural Systems. AIAA/ASME/ASCE/AHS 20th Structures and Materials 
Conference, St. Louis, Missouri, April 1979. 

4. Ormiston, R. A.: Comparisons of Several Methods for Predicting Loads on 

a Hypothetical Helicopter Rotor. Proceedings of the Meeting of Specialists 
on Rotorcraft Dynamics, American Helicopter Society and Ames Research 
Center, February 1974. 

5. Flannelly, W. G. ; Bartlett, F. D., Jr.; and Forsberg, T. D.: Laboratory 

Verification of Force Determination, A Potential Tool for Reliability 
Testing. USAAMRDL TR76-38, January 1977. 

6. Berman, A.; and Flannelly, W. G.: Theory of Incomplete Models. AIAA 

Journal, vol . 9, August 1971, pp. 1481-1487. 

7. Berman, A.: System Identification of a Complex Structure. AIAA Paper 

75-809, Denver, Colorado, May 1975. 

8. Chen, J. C.; and Garba, J. A.: Matrix Perturbation for Analytical Model 

Improvement. AIAA Paper 79-0831, St. Louis, Missouri, April 1979. 

9. Kennedy, C. C.; and Pancu, C. D. P.: Use of Vectors in Vibration Measure- 

ment and Analysis. Journal Aeronautical Sciences, vol. 14, no. 11, 

November 1947. 


145 



10. Klosterman, A. L.; and Lemon, J. R.: Dynamic Analysis via the Building 
Block Approach. Shock and Vibration Bulletin, 42 (4), 97-104, 

January 1972. 

11. Ewins, D. J.: Whys and Wherefores of Modal Testing. SEE Journal, 

September 1979. 


146 



./ 


1. Report No. 

NASA CR-3429 


4. Title and Subtitle 

ANALYTICAL TESTING 


2. Government Accession No. 


3. Recipient’s Catalog No. 


5. Report Date 
May 1981 


6. Performing Organization Code 


7. Author(s) 

W. G. Flannelly, J. A. Fabunmi, and E. J. Nagy 


9. Performing Organization Name and Address 
Kaman Aerospace Corporation 
Old Windsor Road 
Bloomfield, Connecticut 06002 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


8. Performing Organization Report No. 

R-1614 


10. Work Unit No. 


11. Contract or Grant No. 

NAS1-15414 


13. Type of Report and Period Covered 
Contractor Report 
6/27/78 - 9/15/80 


14. Sponsoring Agency Code 


15. Supplementary Notes Langley Technical Monitor: F. D. Bartlett, Jr. 

The Contract research effort which has led to the results in this report was 
financially supported by the Structures Laboratory, AVRADC0M Research 
and Technology Laboratories. Final Report 

16. Abstract 

Analytical methods for combining flight acceleration and strain data with shake test 
mobility data to predict the effects of structural changes on flight vibrations and 
strains are presented. This integration of structural dynamic analysis with flight 
performance is referred to as analytical testing. The objective of this methodology 
is to analytically estimate the results of flight testing contemplated structural 
changes with minimum flying and change trials. The category of changes to the 
aircraft includes mass, stiffness, absorbers, isolators, and active suppressors. 
Examples of applying the analytical testing methodology using flight test and shake 
test data measured on an AH-1G helicopter are included. The techniques and procedures 
for vibration testing and modal analysis are also described. 


17. Key Words (Suggested by Author(s)) 

Structural dynamics 
Flight testing 
Shake testing 
Vibrations 
Absorbers 


19. Security Classif. (of this report) 
Unclassified 


18. Distribution Statement 


Modal analysis 
Stress analysis 
Helicopters 
Suppressors 
Isolators 


20. Security Classif. (of this page) 
Unclassif ied 


Unclassified - Unlimited 


Subject Category 39 



For sale by the National Technical Information Service. Springfield, Virginia 22161 


NASA-Langley, 1981 



















